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FOREWORD 


This report was prepared for the National Aeronautics and Space 
Administration by the Applied Science Associates under Contract NAS6-2307. 
J. T. McGoogan and H. R. Stanley, of NASA Wallops Station, acted as NASA 
coordinators. 

The study was performed both at the Applied Science Associates and at 
North Carolina State University. N. E. Huang served as project director. 
The theoretical part was assisted by Professor C. C. Tung of the North 
Carolina State University at Raleigh as a conslultant , Dr. N. Guttman 
worked on the dispersive relation, and Mr. S. R. Long worked on the 
capillary gravity waves interaction with currents. The experimental part 
was carried out at the North Carolina State University through a sub- 
contract to Professor F. Y. Sorrell. G. V. Strum and S. R. Long also 
worked on the experiments. 
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ABSTRACT 


This report presents the results of a one-year study of current-wave 
interactions. The subjects treated include both theoretical and 
experimental parts. 

In the theoretical study, three problems are discussed. The first is 
the dispersive relation of a random gravity-capillary wave field. It is 
concluded that no universal relationship can be found for all the sea 
states. The second problem is on the changes of the statistical properties 
of surface waves under the influence of currents. The possibility of 
utilizing such changes for remote sensing of surface currents is also 
discussed. The third problem is on the interaction of capillary-gravity 
with the non-uniform currents. 

The experimental study deals with the measurement of wave-current 
interaction and trying to establish the feasibility of using such measurements 
for remote sensing of surface currents. A laser probe was developed to 
measure the surface statistics. In a laboratory, the possibility of using 
current-wave interaction as a means of current measurement was’ demonstrated. 


iii 



CONTENTS 


Page 

FORWARD ii 

ABSTRACT iii 

CONTENTS - . iv 

PART I THEORETICAL STUDIES 

Section 

1 DISPERSIVE RELATION FOR GRAVITY-CAPILLARY WAVES 

1.1 Introduction 2 

1.2 Specification of the random sea 16 

1.3 Mathematical derivation and physics of the dispersive 

realtionship 21 

1.4 Quantitative estimates of the nonlinear dispersion 32 

References 4? 

Appendices 50 

2 STATISTICAL PROPERTIES OF SURFACE WAVES UNDER NON-UNIFORM CURRENTS 

2.1 Introduction 66 

2.2 Distribution of wave elevation 66 

2.3 Distribution of wave amplitude 71 

2.4 Expected number of zero crossings and maxima 76 

2.5 Velocity of zero 78 

2.6 Conclusion 86 

References 90 

3 INTERACTION OF CAPILLARY-GRAVITY WAVES WITH NON-UNIFORM CURRENTS 91 

References . 101 

PART II EXPERIMENTAL STUDIES 

Section 

1 INTRODUCTION 103 

2 WAVE MEASUREMENT SYSTEM 

2.1 Introduction 107 

2.2 Principle of Operation 109 

2.3 The Diode Matrix 112 

2.4 Electronics for the Diode Matrix 113 

2.5 Summary 116 

3 WAVE DECAY STUDY 

3.1 Introduction 129 

3.2. Studies of the effect of surface films 130 

iv 



3.3 Data acquisition techniques 

3.4 Capillary waves on constant currents * 

3.5 Wave decay measurement on a constant current 

3.6 Capillary waves on a variable current 

4 EXPERIMENTAL STUDIES OF CAP ILLARY- GRAVITY WAVE INTERACTION 

4.1 Introduction 

4.2 Experimental technique 

4.3 Experimental results 

4.4 Summary and conclusions 

5 PRELIMINARY WIND-WAVE MEASUREMENTS 

5.1 Introduction 

5.2 Data 
References 


135 

137 

140 


159 

160 
163 
175 


176 


176 

180 



v 



PART I 


THEORETICAL STUDIES ' 

The theoretical study consists of three parts. The first* part 
discusses the non-linear dispersive relation in a random gravity- 
capillary wave field. The result indicates that the dispersive relation 
is a function of the energy spectrum; therefore, a universal dispersive 
relationship seems impossible. The second part deals with the statistical 
properties of surface waves under the influence of currents. The changes 
due to currents can be used as an additional method to measure surface 
current. The third part dicusses the interaction of capillary-gravity 
waves with non-uniform current. A perturbation scheme is used to overcome 
the difficulty of multiple values of the dispersive relationship in this 
range of the wave spectrum. 



1* Dispersive relation for gravity— capillary waves. 


INTRODUCTION 

Looking over the surface of the open ocean, a casual observer notices 
a seemingly endless moving succession of irregular humps and hollows, or 
waves. These surface waves are the most common motion at the air-sea in- 
terface. They form an oscillatory pattern and travel in all directions 
and at different speeds. 

One characteristic of an oscillatory wave is its period. This quan- 
tity is the time interval between passages, at a fixed point, of a given 
phase of the oscillation. Its reciprocal is the wave frequency n. Since 
n is a function solely of time, it can be measured relatively easily from 
a point by noting the elapsed time between the passage of two crests. 

Another characteristic is the wave length, or the distance between 
particles moving in the same phase of the oscillation. The reciprocal 
of this quantity is the wave number k. Since the wave travels on the two 
dimensional free surface x = (x,y), the wave length and directional prop- 
erties can be combined and represented by the vector wave number 
IS = Measurement of the directional wave length at infinitely many 

points is necessary to obtain spatial correlations over the two dimen- 
sional wave number plane. 

The velocity c at which the wave travels is defined as 
„ _ n 

? = X 0) 

Since the frequency n is non-di rectional and relatively easy to measure, 
and since the wave number k is economically and logistically prohibitive 
to obtain over the free surface, it would be advantageous to find a re- 
lationship between n and k. 
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For a single gravity wave traveling in the x direction, such a re- 
lationship can be derived. Consider a two dimensional spatial coordinate 
system (x,z) where z is measured positively upward and the plane z = 0 
represents the mean free surface level as shown in Figure 1. Assuming 
an incompressible, inviscid fluid, the motion can be regarded as irrota- 
tional and is thus controlled by Laplace's equation 

V <f> (x,z,t) « 0 (2) 

where $ is the velocity potential and t is time. By requiring that no 
motion exist at z = -h, where h is the water depth, and that Bernoulli's 
linearized equation 


9 ? 


3t p 


(3) 


hold at the free surface, (2) can be solved for 1 In (3), g is the 
gravitational acceleration, p is the pressure at the interface, p is the 
density of the water, and c is the oscillatory surface elevation 


c(x,t) = a cos(kx - nt) 


(4) 


where a is the wave amplitude, Bernoulli's equation has been linearized 
under the assumption that terms of the order of the square of the veloc- 
ity components are negligible compared to other terms. 

The solution to (2) is 



cosh[k(z+h)] 

sinh(kh) 


sin(kx - nt) . 


(5) 


Substituting (5) into (3), using (1), and assuming that ke«l, , a 
small wave amplitude assumption, 
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[g - c 2 k coth{kh)]c = - \ | . . ' (6) 

If capillary waves are also included, the effect of- surface tension 
must be considered. The difference between the water pressure p and the 
atmospheric pressure p a at the interface is a result of surface tension 
such that 

2 

p = p - PY 14 (7) 

3X 

where y is the ratio of the surface tension to the water density. Equa- 
tion (7) has been 1 inearized by considering terms of the order of the 
square of the gradient of the surface elevation in the representation of 
the surface curvature to be negligible compared to other terms. The at- 
mospheric pressure is assumed to be constant, and the free surface c * 0 

is taken as a reference level at which p * 0. In this situation p. = 0 

a 

since the surface tension term vanishes, Thus, 

2 

5 C 

P B • PY — z • (8) 

ax 

Substituting (4) into (8) and then (8) into (6) and rearranging, 

2 2 

g + yk “ c k coth(kh) . ( 9 ) 

For deep water waves, h is large and asymptotically 

coth(kh) *= 1 

so that 

£ + Yk ■= c • - (11) 

Equation (11) is the linearized relationship for the speed of a 
single wave under the restoring Influences of gravity and surface tension. 
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It is depicted graphically in Figure 2. The speed is a function of the 
wave length and thus represents the dispersion of the wave. Substituting 
(1) into (11) results in the dispersive relationship as a function of 
n and k. 



3 

gk + yk 


(12) 


Equation (12) is graphically depicted in Figure 3. 

It is important to note that (12) is in reality only approximately 
true since just linear terms are considered. Cole (1968) discusses the 
general case of periodic motions of slightly non-linear oscillations . 

He treats various physical problems that are characterized by the pres- 
ence of a small disturbance which acts over a long time. He shows in 
his non-linear asymptotic expansion of the solution to these oscillatory 
problems that a frequency shift from the linear solution will occur, and 
that this shift is generally a function of the amplitude of the oscil- 
lation. 

Becoming more specific. Kinsman (1965) discusses the non-linear ex- 
pansions of the phase speed for single gravity and capillary waves. For 
gravity waves, the third order expansion includes a perturbation on the 
linear phase speed that is a function of the square of the amplitude. 

The capillary wave higher order solution also includes a perturbation on 
the linear phase speed that is a function of the square of the amplitude. 
Using (1) to transform these results into frequencies, it is found that 
even though the restoring forces to the oscillatory waves are different, 
both the gravity and capillary wave frequencies will exhibit shifts from 
the linear solution that are dependent upon the square of the amplitude 


of the waves. 
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Comparison of the linear dispersion with observed data was made by 
Longuet-Higgins et al . (1963) for gravity waves with frequencies up to 
4 rad sec . They found that the waves appeared a little shorter than 
those predicted by linear theory. The mean increase in the observed wave 
number is about 10 percent. The authors note, however, that they are 
skeptical of the calibration of the instruments. In addition, the buoy 
used has a diameter that is about one-third the shortest wave length. 

When the buoy is located near the crest or trough of the shorter waves 
in the range of the experiment, it will not completely rest on the water 
surface and slope measurements may be inaccurate. The resulting wave 
numbers may therefore be inaccurate. These possible errors cast some 
doubt on the validity of the 10 percent difference between the observed 
wave numbers and those suggested by linear theory. 

Yefimov and Khristoforov (1971) also compared the linear dispersion 

with observed data for gravity waves with frequencies up to 4 rad sec“\ 

Their comparison is based on measurements of the velocity spectra at 1 m 

and 3 m below the surface. They show that for frequencies below 2.5 
-1 

rad sec the waves are slightly shorter than those suggested by linear 
theory. Above 2.5 rad sec , the waves are longer than suggested by 
linear theory. They emphasize, however, that small deviations from the 
linear relationship in the low frequency range can be caused by the re- 
cording by the measuring equipment of harmonics of the principal oscil- 
lation. Additionally, the authors point out that the large deviations 
between the observed data and the linear theory at frequencies near 4 
rad sec ^ are the result of the high frequency turbulent regime of the 
velocity spectra. 
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Grose (1971) and Grose et al. (1972) measured wave heights relative 
to a buoy for gravity waves with frequencies up to about 4 rad sec”^. 

It is shown in these papers that the linear dispersion suggests wave 
numbers that are too large when compared to values computed from the 
wave height measurements. In other words, the observed wave numbers 
are smaller than those predicted by linear theory. The authors state, 
however, that the measuring technique does not differentiate between 
fundamental waves and harmonics. All wave components are treated as 
fundamentals so that errors in the observed dispersive relationship do 
exist. 


Comparison between the linear and observed dispersion in the capil 
lary wave range has, unfortunately » not been made. Experimentation in 
this area needs to be accomplished AnH chrml rl nv'nwiylA vov'vr >•!« 1 a 

| " ^ ^ " * ** < W | 1 K* t f II W I L# I J VY I J | | C 


results. 

Linearizing the governing equations and boundary conditions de- 
creases the reality of the surface motion but greatly simplifies the 


mathematics since the principle of superposition is applicable. Any 
linear combination of solutions to the equations is itself a solution. 

For example, consider two waves of equal amplitude but slightly differ- 
ing frequencies and wave numbers traveling in the same direction. Mathe- 
matically, 


a a 


cos(k x - 

l 


n t) 
i 


(13) 


C = a cos(k x - n t) 
2 2 2 


( 14 ) 



05 ) 


11 


k - k , n ~ n 
1 2 1 2 


k - k « k , n-n «n 
12 112 1 


06 ) 


Each wave travels at its own phase speed. When they meet, however, they 
interact to form a new wave. 


C c C + C 

3 12 


= 2a cos 


(k -k^x (n - n^t 


cos 


(k+k 2 )x 


L 


(n 1 +n 2 )t 


07 ) 


From (15) and (16), the argument of the first cosine function of (17) will 
be small. This term thus serves as a slowly varying amplitude modulator 
of the second cosine function or basic wave. Writing 

(k -kjjx (n 2 -n 1 )t 


and 


a = 2a cos 


k 1+ k 2 

k 3 = g * n 3 


n i +n 2 


08 ) 


(19) 


(17) can be restated as 


r = a'costk.x - n t) . 


( 20 ) 


Since k 3 and n 3 must satisfy the dispersion depicted in Figure 3, 


2 3 

n * 9k 3 + Yk 3 


( 21 ) 
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or 


n +n 2 k +k k +k 3 

(-i-i) . g + yf-lj-i) . 


( 22 ) 


Obviously, from (22), even though the new wave was produced from the sum 
of and c 2 » ^ n °t be dispersed according to the summation rela- 
tionship 

22 , 3 3 

n ! + n 2 ■ gCkj+l^) + Y(kj+k 2 ) . (23) 


The addition of the two sinusoids as seen in Figure 4 begins to 
show some of the properties of ocean waves. By adding together a large 
number of sinusoids and by considering the directional properties of the 
waves over the two dimensional sea surface plane, the irregular pattern 
of the ocean surface can be approximated. The true random surface, how- 
ever, must he described by considering both the linear and the non-linear 
parts of the governing equations. 


The importance of non-linear interactions for gravity waves was dis- 
cussed by Phillips (1960a, 1969). He showed that the interaction of 
three primary waves can cause a resonance condition to exist so that 
secondary waves will be generated. A weak energy transfer from the pri- 
mary to the secondary waves will also occur, McGoldrick (1965) studied 
interactions among capillary and among capillary and gravity waves and 
found that resonance and thus secondary wave generation and correspond- 
ing energy transfer occur with certain triads of waves. Longuet-Higgins 
(1963a) explained the non-linear mechanism for the formation of capil- 
lary waves on the forward face of gravity waves. He related this process 
to the generation of waves by wind. 
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All of the above studies are important contributions to the under- 
standing of surface waves. They all concern non-linear Interactions, but 
the development of each theory relies on the assumption of the validity 
of the linear dispersive relationship (12). The resonance models in 
particular are critically dependent upon this assumption. The reliance 
upon the linear dispersion is true not only for the three papers cited, 
but for virtually all studies involving non-linear interactions. If (12) 
Is not a good representation of the true dispersion, then all of these 
papers — some of which are major contributions--may be improved. 

The present study investigates the dispersion of non-linearly in- 
teracting random pelagic gravity and capillary surface waves. A quanti- 
tative analysis is given, and the deviation from linear theory is exam- 
ined in terms of the non-linear interactions. 

Not only is the non-linear dispersion important in the assessment 
of major theories that lead to the understanding of random waves, but it 
is also important in determining the transformations between wave number 
and frequency spectra. An accurate transformation is necessary in view 
of the economical and logistical problems of wave number measurement in 
the random ocean and of the relative ease of frequency measurements. 

The discussion so far has related to wave-wave interactions. Wave- 
current interactions, however, also affect the dispersion. The change in 
wave number and amplitude when short waves ride on long swells was rigor- 
ously calculated by Longuet-Higgins and Stewart (I960*). The changes in 
wave characteristics when waves interact with currents was also noted by 
Whitham (1960, 1962). Since an accurate transformation between wave number 
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and frequency spectra is desired, the present study additionally describes 
the modification needed in the dispersion when the waves are under the in- 
fluence of current. Although the modification is given only for a single 
wave, the logic of the derivation is valid for random waves. 
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1.2 SPECIFICATION OF THE RANDOM SEA 

The detailed pattern of the waves on the sea surface is random in 
the sense that it varies without regularity in both space and time. The 
motion of the waves obeys certain known dynamical and kinematical restric- 
tions. The randomness, however, must be treated statistically. The sta- 
tistical measures of the motion can then be used to ascribe some order to 
the surface structure. It is necessary, therefore, to specify the dynam- 
ical, kinematical, and statistical properties of the random sea surface. 

Consider a rectangular coordinate system with the z-axis measured 
positive vertically upward with the plane z = 0 at the mean free surface 
level. Any point in this system can be described by z and a horizontal 
position vector x = (x,y). Under the standard assumptions of an incom- 
pressible, in viscid fluid, the motion can be approximated as irrotational . 
This motion is governed by Laplace's equation 

2 

V *(x,z,t) = 0 ( 24 ) 

where <f>(x,z,t) is the velocity potential as a function of x,z and time t. 

If it is further assumed that the random wave field is statistically 
stationary with respect to both time and space and that there is no 
motion at z - Phillips (1960a) showed that the solution to (24) is 

4 >(x,z,t) = ™ dA(k,n)exp[|k| z]exp[i(k-x - nt)] ( 25 ) 

where dA(k,n) is any complex valued random function of the horizontal 
wave number vector k and frequency n. The velocity potential is rep- 
resented as a Fourier-Stiel tjes integral with the integration over all 
wave number-frequency space. 
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Under the same assumptions Phillips (1960a) showed that the surface 
elevation c(x,t) can be represented as 

c(x,t) = j^dB(k,n)exp[i(k*x - nt)] (26) 

where c is the displacement from the mean free surface c* = 0 and dB(k,n) 
is another complex valued random function. Since s(x,t) is real, 

c(x.t) = j^dB(k,n)exp[i(k*x - nt)] 
e t * n ) ex PC"i ( is - nt)] 

= ^dB*(-k,-n)exp[i(k«x - nt)] (27) 

where dB*(k,n) is the complex conjugate of dB(k,n)* t 

A relationship exists between dB(k,n)- and the directional wave energy 
spectrum X(k,n) such that 


dBfkj )dB*(kj ,nj) 


X(k. ,n^ Jdk^ dn* 
|^0 otherwise 


if 1 = j 


(28) 


where the overbar indicates an ensemble average. The spectral function 
Is perhaps the most common representation of the wave field. In order 
to relate quantities other than the surface elevation to the wave spec- 
trum, it is necessary to find a relationship between dA(k,n) and dB(k,n). 
This can be done from the kinematic boundary condition at the free sur- 
face . 
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If the position of the free surface is specified by e» the total 
derivative of c is 


ck_ = . 

dt 3t 


vv 


( 29 ) 


where = (a/ax, a/ay) is the horizontal gradient operator, q is the 
horizontal velocity vector, and the subscript c indicates measurement 
at the free surface. For irrotational flow (29) becomes 


fi = (8*/az) - (v h *) c -v h ; . 


(30) 


Equation (30) represents the free surface kinematic boundary condition* 
Substituting (25) and (26) into (30), 

j^-ind3(k,r.)exp[i (k*x -nt)] - |kj uA{k,n)exp[ jkjc]exp[1{k*x - nt)] 

+ ^ k-k 1 dA(k»n)dB(k 1 ,n 1 ) exp£ | k | c]exp{i [( k+k a ) -x - (n + n 1 )t]} 

( 31 ) 

The exponential term in (31) involving the surface slope |k|c can be 
expanded in a Taylor series around c = 0 such that 

i 2 i • 


exp[|k|c] = 1 + | k | c + |k | c] + ... yjC 1 k I + ... + R. (32) 


The convergence of the series for gravity waves can be determined by 
noting that the maximum steepness a gravity wave can have before it 
breaks is .142. This value was determined by Stokes (1880), and accord- 
ing to Kinsman (1965), it may be taken as an established value. In terms 
of wave slope, the maximum value is one half the maximum steepness, or 
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.071. Terminating (32) after only two terms results in a maximum error 
R 


1 2 

R - £{.071) exp( .071 ) * .003 . (33) 

This error is negligible. 

For capillary waves, Crapper (1957) found that the maximum steep- 
ness can be .730. The maximum slope is therefore .365. The maximum 
error if only the first two terms of (32) are used is then 


1 2 

R = £<.365) exp{ .365) = .096 . 
If the first three terms are used, 


(34) 


R = (l/6)( .365) exp( .365) - .012 . 


(35) 


Although these errors are larger than that for 


ytvmuj 


na vco , Li 


hey ai 


4 C 


still small. 

Using the series expansion for exp[jkjs] truncated after the first 
three terms, the relationship between dA(k,n) and dB(k,n) is found by 
successive approximations. Details of the operation, following the 
method developed by Huang (1971) in his second order study of Stokes 
drift in a deep water random gravity wave field, are in Appendix A. To 
the third order the relationship is 


dA(k,n) 


- jfcy ^B(k ,n) 
4B(kj + i 


+ i 


«/ / 


!h n ! 


|k||k - k,l 


/ / / / (n ' n I - n 2 } 

I *1 


n - n 1 )dB( k - k ,n - n ) 
2^- + -l _ S 2 ) 
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-k-(k - k 1 ) (k - kj).(k - ^ - k 2 ) 
|k - k,||k - k - k | 


dB(k - kj - k 2 , n-rij - n 2 )dB(k 1 ,n 1 )dB(k 2 , n 2 ) 


A dynamic boundary condition on the free surface is imposed by the 
requirement that the pressure on the two sides of the surface differ only 
as a result of surface tension. From Bernoulli's non-linear equation, 
the pressure in the water at the free surface is given by 

n l 2 • 

f+ 9C + (a*/8t) c + = 0 . (37) 

The difference in pressure resulting from the surface tension can be 
written as 

P = P a - py(v^c)[l + (v h c) 2 J-< 3 / 2 > (38) 

where the second term on the right represents the contribution from the 
surface tension. As in the linearized equations, the atmospheric pres- 
sure is assumed to be constant. If the free surface c = 0 is taken as 
a reference level at which p = 0, then p a = 0 since the surface tension 
component vanishes. The resulting dynamical free surface boundary con- 
straints on the wave motion are 

2 2 2 

gc = -(a*/at) c - g{v*) t + y(v h c)[l + (7 h c) ]'( 3/2 ) . (39) 


t - r 
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1.3 MATHEMATICAL DERIVATION AND PHYSICS OF THE 
DISPERSIVE RELATIONSHIP 

It has been stated that the dispersive relationship is one particular 
property of gravity and capillary waves relating wave number and frequency 
and that several investigators use the relationship in a form derived from 
linear theory for a single wave. The present analysis describes the mean 
deviation and the scattering from the simple dispersion for the more real- 
istic random, non-linearly interacting waves. The analytical derivation 
makes use of the full non-1 inear equations for a capillary and gravity 
surface random wave field and is an extension of the study by Huang (1972) 
on gravity waves. 

The free surface dynamical boundary constraints yield an identity 
from which the dispersion for a single wave train was found.. They also 
provide a corresponding relation for a random wave field. Looking at 
(39), the term raised to the -3/2 power can be expanded into the binomial 
series 

[1 + (V) 2 ] 3/2 “ 1 - ( 3 /2)(V) 2 + 05/8)(v h c) % + ... 

■. [-(3/2) - j + l][-(3/2) - J + 2] ... r-(3/2) - 11 (- 3 / 2 ) 

j! 

(v h c) J + ... + R . (40) 

Equation (40) is valid for 


i<v )2 i < 1 • 


(41) 
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This condition is always satisfied in light of the maximum wave slopes 
permitted before breaking occurs. Truncating (40) after the first two 
terms yields a maximum error 

R = 05/8)(v h c) 1 * . (42) 

Since the maximum gradient is small in order to retain the wave stability, 
the maximum error after a two term truncation is minute. 

Substituting the first two terms of the binomial series (40) and the 
first three terms of the Taylor series (32) in (39) and then using (25) 
and (26) to substitute for $ and c, an expression is obtained in terms 
of dA(k,n) and dB(k,n) such that 

^gdBfk.nJexpfi (k-x-nt)] = j^1n[l + | k | c + j J k J c ]dA(k,n)exp[i (k*x - nt)] 

" \ knk^^J-l l!Sx I” + dk| + lk 2 1 )c] 

dA(k,n)dA(k 1 ,n i )exp{i[(k + - (n + n 2 ) t]} 

rr . 2 

“ T kn I ~ i dB (M)exp[i(k*x " nt )3 ~ y(3/2) 

2 

tak n k { HSI *Ji*!S 2 dB (!S»n)dB(k 1 ,n 1 )dB(k 2 .n 2 ) 

~ 1~2 2 z z 

- exp{i[(k + k + k )*x - (n + n, - njt]} 

~ ~ 1 -Z 12 


( 43 ) 



23 


Substitution of (36) in (43) results in the non-linear third order dis 
persive relationship for random gravity and capillary waves 

S + y|(s| 2 - Tjrr = £ „ f^Bfk^nJexptifk-x - n t)] 

1 - 1 -l l 1 


‘ k,£ k n f 2 dB(ki,n 1 )dB(k 2 ,n 2 )exp{iC(k 1 + k 2 ) 


•x-(n 1+ n 2 )tJ} - 2-k,n,k n - k n - n ) 

-11-22 * 


dB {k 2 .n 2 )e xp [i(ki-x - rijt)] 


(44) 


where 


f = 


2 

n - 


k*(!$ + !<,) i 

i , , i. ■ .i ■ .i i . i. ri ( n -f « ^ a JL /l 

| k j j k + k 'V ' 2 v ' 


k-kj 


Tkn^r /,,n 


/ 

V*+3; 


f. * 


n(n + n + n ) 


1 2 


2 ~ ] k + + k 


2' L‘ 


1, , (k + k + k ) *(k + k ) (k + k ) -k 

& * ^ li. yit i ~ ' 


(k + k )-k 

+ ~W\ — 

(l!$l + IkiDwij + J 


1 o 1 F k “k 

n(n + V - jn |k| - [ [l - 


(k + 1 (!< + k 2 )* k 

1 ' |k + k 2 ||ki| WT^ nn i 


* (3/2)y ! k I k^k,. 


f 3 = 


1 - 


k-k x 


TYTTI7 


~ n(n - n ) 
Ikj - k 2 l 


(46) 

(47) 


Details of the mathematics of the derivation are given in Appendix B. 
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Equation (44) is the most general form of the dispersive relation- 
ship for a random sea surface through third order interactions. It is a 
function of both space and time. If all the non-linear terms on the 
right hand side are neglected, and if only a single wave train is con- 
sidered, the frequently used simple relationship (12) is recovered. 

The deviation from the linear theory is given by the right hand side 
of (44). The mean deviation is easily computed by taking the mean of 
each term of the deviation. The first term represents a random scatter- 
ing of the surface elevation and has a zero mean. 

It is seen from (27) that 

dB(k ,n) * dB*(-k ,-n) * (48) 

The mean Fourier-Stieltjes components In the third term can then be 
written as 


dB(k - k ,n - n )dB(k ,n ) = dB(k - k ,n - n)dB*(-k ,-n ) 

-i -2 1 2 '2 2 -1 -2 1 2 "2 2 


(49) 


and from (28) 


dk^ - k 2 7n 1 “~n 2 )dB*(-k 2 V-n 2 ) = X(k £ ,n 2 )dk 2 dn 2 if ^ = 0 


0 otherwise. 


(50) 
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If k x * rij = 0, then (47) becomes zero. The mean of this third term 
therefore vanishes. 

Equations (48) and (28) can be used to determine the mean of the 
second term. Thus, 


dB ( k x ,n 1 )dB*(-k 2 ,-n 2 ) = A(kj * n 1 )dk 1 dn 1 if ^ + k 2 = n : + n 2 = 0 


0 otherwise . 


(51) 


Substituting k = -k 2 , = -n 2 into (46) yields a non-zero term 

(k + k>k 


f 2 = ]kT 


,, k-(k - k.)(k - k.)-k 

2^~ + 2 -P‘k |k - kj [Ff 


. ^ T n (n + n.) 

“Til 


|n 2 |k| - j 


1 - 


jc-kj 




(|k| + |k 1 |)nn 1 


♦i 


1 - 


(k - kp-kj “ 

(k - kp.k 

1 

1 

t 

JSt 1 

1 

|k| 


I !S 1 2 MS i 1 2 • (B2) 


It is seen that the mean of the second term is therefore non-zero. 
The resulting mean non-linear dispersion is 

2 n 2 . / / 


9 + v|k| - yjjj- = -k i n i f 2 X( !? 1 »n 1 ) d !s l dn i 


(53) 


The right hand side of (53) represents the mean non-linear deviation 
from linear theory in terms of the spectral function. This expression 
is greatly simplified if only unidirectional waves are considered. 
Equation (23) reduces to 
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2 ^ 

9 + Yk - -f- = ■-{ { C 3 n 2 k l + nn^k + kj) - (3/2) Y k 2 k 2 ] 

f X(k x ,n 1 )dk 1 dn 1 . (54) 

Since the spectral function is even and n^, are odd functions, (54) 
becomes 

2 

9 + yk ' *F* “ 'k n t nn i k i - (3/2) Y k 2 kj]X(k 1 ,n 1 )dk 1 dn 1 . (55) 

The mean squared random error is a measure of the variability of 
the random scattering about the mean deviation. The right hand side of 
(44) represents the deviation e from the linear theory. The mean of its 
square , is the mean squared random errors The computation of to 
the third order requires that we consider only the first term on the 
right hand side of (44) as inclusion of the other terms will yield high- 
er than third order results. Thus, 

£ S 'k 1 n 1 f i dB ^i* n i^ exp ^ 1 * " n i t ^ * ■ (56) 

The mean squared random error is arrived at by once again using (48) 
and (28) to compute the mean of the square of (56). The resulting ex- 
pression is 

T 2 " • ( 57 ) 

For unidirectional waves, (57) reduces to 
p = n 2 k 1 n 1 n > X < l V n i> dk i d V 


( 58 ) 



Anticipating the results of the next section, the behavior of the 
mean deviation and root mean squared random error can be schematically 
visualized as in Figure 5. The ratio of the linear dispersion to the 
non-linear dispersion is plotted against the wave frequency. It is seen 
that as n increases from the low frequency gravity waves through the 
high frequency capillary waves, the mean deviation and random scatter- 
ing first increase in magnitude and then decrease in magnitude. The in- 
flection points of the curves are found in the transition zone between 
the gravity and capillary wave ranges. 

The reversal of the mean effects can be seen by looking at the term 
in brackets in (55). The first component dominates in the gravity wave 
range since in this range the square of the wave number is much less than 
the wave number itself. In the capillary wave range, however, the wave 
number is large so that its square is even larger. The second term 
therefore dominates in this range. Since the two terms are of opposite 
sign, their effects will be opposite. 

The random fluctuations are a function of the quasi -Euleri an veloc- 
ity components. The schematic depiction shows that these random fluc- 
tuations are maximized in the transition zone between the gravity and 
capillary waves. This curve is based in part on measurements of veloc- 
ity spectra by Yefimov and Khristoforov (op.cit.). These authors show 
that in the gravity wave range the random, turbulent type velocity com- 
ponents become increasingly stronger as the transition zone is approached 
No measurements of velocity spectra have been made in the capillary wave 
range, but Phillips (1969) discusses integrated frequency spectra in this 



gk + 7 k 3 


GRAVITY WAVES 


TRANSITION ZONE 


CAPILLARY WAVES 


Figure 5. Schematic of the effect on the linear dispersion of the non-linear mean deviation ( 
squared random error ( *— — - - ) 
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range. He states that the maximum spectral values occur in the lowest 
frequencies, or those frequencies near the transition zone. At higher 
frequencies, the spectral values decrease rapidly because of the in- 
creased effect of viscous dissipation of the capillary waves. Thus, as 
the frequency increases, the wave energy becomes less and therefore the 
velocity components decrease in magnitude. 

The dispersion by definition relates the velocity at which a wave 
travels to the wave length. This has the effect of sorting the waves 
as they travel. For simplicity, consider three gravity waves of dif- 
fering lengths propagating unidi rectionally from a point x Q at time t . 
Since the waves are traveling at different speeds, at a later time t 1 
the fastest wave will at at x 1 , the next fastest at x 2 <x : , and the slow- 
est at X 3 <X 2 <X 1 * It 1$ therefore seen that even though the three waves 
started traveling from the same point at the same time, they will be 
sorted at t r 

The linear dispersion (11) shows that the square of the phase speed 
Is indirectly proportional to the wave number when the wave is restored 
to the mean free surface 'level by gravity and directly proportional to 
the wave number when the wave is restored by surface tension. The ef- 
fect of non-linear wave-wave interactions can be readily seen by look- 
ing at the simplified case of unidirectional waves. Using (1), the mean 
non-linear dispersion (55) can be rewritten as 

cZ = -f + yk + n ) k i - (3/2) Y kk*]X(k 1 .n 1 )dk i dn 1 
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* c o 


1 + if kiVi^i-V^V (3/2) 2 

11 C 

0 u o 


k 1 n 1 k i X(k i* n i )dk i dn i 


(59) 


where c 2 represents the linear dispersion. If interactions did not oc- 
cur, c 2 would be coincident with c^ and the linear theory would suffice. 
Physically, then, the third order mean non-linearities behave in the 
sense of a deviation from the linear dispersion. They thus affect the 
dispersion by changing the phase speed c Q that results from linear theory, 
Looking at the mean deviation, the first integral on the right hand 
side of (59) is an expression for the mean random ocean quasi-Eulerian 
velocity as derived by Phillips (1960b). This velocity is the mean of 
the horizontal component of the fluid velocity at the free surface for 
a fixed point (x,y). It is a quasi-Eulerian quantity in that it is meas- 
ured at a fixed point on the projection of the free surface, but the 
vertical coordinate is allowed to move up and down with the free sur- 
face. A similar formula was derived by Longuet-Hi ggins and Phillips 
(1962) from a generalization of the interactions between two waves. The 
difference between the two expressions is a factor of two. Phillips 
(1960b), however, showed that the mean Eulerian velocity is one-half the 
Stokes drift. The Stokes drift is the mean velocity following a fluid 
particle and is by definition a lagrangian property. Since the present 
study is based on an Eulerian reference frame, the first integral appears 


reasonable. 
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The second integral in (59) represents the mean square of the ratio 
of the surface elevation to the wave length. It is multiplied by the 
surface tension and is thus important in the capillary wave range. This 
term also appears reasonable in view of Crapper's (op.cit.) exact solu- 
tlon to the finite amplitude non-linear single capillary wave problem. 

He showed that the square of the phase speed is influenced, as in the 
present study, by a term proportional to the square of the ratio of the 
surface elevation to the wave length. 



1.4 QUANTITATIVE ESTIMATES OF THE NON-LINEAR DISPERSION 


The effect of the mean non-linear deviation can be quantitatively 
assessed if the spectral function is known. An order of magnitude esti 
mate can be obtained by using the equilibrium range spectrum. For grav 
ity waves assume the equilibrium spectrum given by Phillips (1958), 

X(k,n) = 6(k - (60) 

rr u •• 

for 

= (g/v) 1/2 . n«n^ = (4g 3 / Y ) 1/4 (61) 

where e is a universal constant and i is the Dirac delta function. The 
constant 3 has been measured experimentally by Pierson (1960), Burling 
(1959) and Longuet-Hi ggins et. al . (op.cit.). An average value is 1.2 x 
10 . For capillary waves Phillips (1969) gave an equilibrium spectrum 

based on dimensional considerations such that 

. X(k,n) = - B - y - 7 y 3 5 (k 3 - ) (62) 

n ' T 

for 


k <<k<<k , n <<n<<n (63) 

where b" is another constant, and k^ and are the cut-off wave number 
and frequency, respectively, at which viscous effects dominate. Cox 
(1958) measured a constant for frequency spectra of the wave slope. 
Phillips (1969) showed that this constant could be used to approximate 
B' as 10“ 2 . Wu (1972), however, showed that this value may be in error 
by a factor of two. 
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In the gravity wave range the surface tension is unimportant com- 
pared to gravity so that after substituting (60) into (59) 



c 2 [i 

°L c 2 


7 7 n*“^k dk dn 

kjn^i V K l n i 


(64) 


where 

n 2 = gk x . (65) 

Using (65) and integrating (64) over the gravity wave range and noting 
that 


c 2 = g/k (66) 

then 

c 2 a + ...eg 9 -" 

C ~0 2 
c n 

L o i A 

= c 2 [1 + 3 n/n L ] . (67) 

Equation (67) shows that the mean phase speed of interacting waves will 
be less than that predicted by linear theory, and that the difference 
will depend on the ratio n/n^ 

Examining this result in terms of wave number and frequency, sub- 
stitute (1) and (66) into (67) and rearrange so that 

gk = n 2 [l + (an/nj)] -1 . (68) 

Expanding the multiplier of n 2 in a truncated binomial series, 

[1 + (Bn/n^]** 1 ^ 1 - Bn/n 1 . 


(69) 
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Hence 

gk - n 2 (l - Bjj—) . 


( 70 ) 


The error in (69) is negligible since n and n are of the same order and 

* 

3 is small. Equation (70) clearly shows that for interacting gravity 
waves, an increase in wave number results from an increase in frequency 
below that predicted by linear theory. Furthermore, the effect of the 
interactions increases as the ratio n/n x increases. This result agrees 
with Huang (1972) and is consistent with Stokes' (1847) theoretical 
solution for gravity waves. Additionally, the result agrees with the 
observed data discussed earlier. 

In the capillary wave range the gravity terms are unimportant com- 
pared to the surface tension terms so that after substituting (62) into 


( 59 ), 


C 2 = c 2 {l + (B7c 2 )^ i (c Y 2/3 n l * 4/3 k 1 - (3/2) Y S/S kn 1 " 7/3 k^] 


dk 1 dn 1 > 


(7V) 


where 


2 ,3 


n. = 




(72) 


Using (72) and integrating (71) over the capillary wave range, 

c 2 = c 2 {1 - (b'/c 2 ) [(3/2) Y kln - 3c Y 1/3 n 1 1/S ] > . (73) 

The mean effect of the interactions on the phase speed can be assessed 
after the sign of 
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(3/2)ykln rij - Zc^^n^* 


(74) 


Is determined. Using (1), dividing by 3y and rearranging (74) is re- 
written as 




(75) 


The capillary wave range in which the equilibrium spectrum is valid 
begins at n^ = 100 sec" 1 as shown by Phillips (1969) and at « 3.6 cm’ 1 
as shown by Wu (op.cit.). . The upper limit to the range is determined by 
the energy dissipation resulting from molecular viscosity. Wu (op.cit.) 
showed that the cut-off wave nuinber k ' :: t' 13 c|" , ; The corresponding 
cut-off frequency n^ - 1,000 sec -1 . Additionally, is the lower bound 
of the equilibrium range so that n x is always less than or equal to n. 
using these limiting values, the first term in (75) varies from .50 to 
.35 as n a increases and the second term varies from .44 to .34 as n in- 
creases. Equation (75) and therefore (74) is thus positive. It should 
be noted, however, that as the viscous dissipation of the wave energy 
becomes more pronounced, i.e. f at the upper limits of the range for which 
the equilibrium spectrum is valid, the difference between the two terms 
becomes smaller and smaller. Since (74) is positive, the square of the 
phase speed c in (73) will be less than that for the situation in which 
no interactions are involved, but only by a small amount. 

In terms of wave number and frequency, for capillary waves 


c 


yk . 


(76) 



Substituting (1) and (76) into (73) and rearranging. 


1 / 3 , 


- (3/2)s'ln n 1 3 + Y *' ''k{3e'nn 1 1/3 ) - n = 0 


( 77 ) 


1/3 

This equation is a cubic expression of y k and can be solved in terms 
of n by using the formulas and procedures analogous to those outlined in 
Appendix C. Cubing the solution results in 

.3 o x or/ '1/2.2/3 l/ 21/3 1/2.1/3. 1/2.2 / 3 

yk s 2r + 3[(r + D ) (r - D ) + (r + D ) (r - D ) ] 


(78) 

where 

r * j n 2 [l - (3/2)e'ln n^" 1 (79) 

and 

D = e'Vn^l - (3/2)e'ln n^" 3 + jo^l - (3/2)e'ln 

Since e" is such a small quantity, the first term in (80) is 
compared to the second term. Thus, D is approximated by 

D * ^“[1 - (3/2)p'ln r^]' 2 . 

Rewriting (78) and using the approximated value for D, 

Yk 3 * 2r + 3(r 2 - D) j/3 [(r + 0 1/2 ) 1/3 + ( r . D 1/2 ) 1/3 ] 

*= n 2 0 - (3/2)e'ln r^]" 1 - 3e'nn 1 1/3 [1 - (3/2)B'ln n 1 ] _1 


2 • (80) 
negligible 


(81) 
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« n 2 n - (3/2}B'ln i^]" 1 - - (3/2)6'7n n l ]’ 4/3 

= n 2 {[l - (3/2)e'1n Dj]' 1 - 3e'(n 1 /n) 1/3 [l - (3/2)e'ln n 1 ]' 4/3 

(82) 

Expanding the terms in brackets into binomial series which are truncated 
so that only terms of the order of or larger are retained, 

yk 3 = n 2 [l + (3e72)ln - Sp^nj/n )^ 3 ] . (83) 


The error from the truncation of the binomial series is of the order of 
2 

. This is negligible compared to the other terms in (83). It has al- 
ready been ascertained that 


1 1/3 

^ In n 1 - (n x /n) > 0 


(84) 


since in the capillary wave range the logarithmic. part will always be 
greater than unity while the ratio will always be less than or equal to 
unity. Therefore, from (83) it is seen that if the frequency is increased, 
the non-linear interactions will have the mean effect of increasing the 
wave number above that which is predicted by linear theory. Furthermore, 
the effect becomes greater as n t and the ratio n/n^ increase. 

The mean squared random scattering defined by (58) can be inter- 
preted as a measure of the random deviation about the mean scattering. 

It is analogous to a variance in the statistical sense. Its square root 
is analogous to a standard deviation. A large value indicates a wide 
spread of data points about the mean while a small value indicates a 
close cluster of points about the mean. 
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Equation (58) can be investigated by looking at the integral func- 
tion. This term represents the mean squared surface velocity spectrum, 
where’the velocity is equivalent to the quasi-Eulerian velocity discussed 

earlier in conjunction with the mean scattering. The root mean squared 
, 1/2 . • 

error, (e 2 ) , then represents a mean of the individually unpredictable 

random variations of the quasi-Eulerian velocity data points about the 
mean quasi-Eulerian velocity multiplied by the constant n. 

The unidirectional non-linear effect on the square of the phase 
speed can be represented in terms of the mean deviation from linear 
theory plus the random scattering. Combining the square root of (58) 
with (55) and then substituting (1) into the combination, 

c 2 e c 2 {1 + (c/c 2 )£ jj k 1 n 1 X(k 1 •n 1 )dk 1 dn 1 - (3/2) 1 -- 

c o 

i n 2 X(k ,n )dk,dn 
ii 

Quantitatively, all components of (85) except the last one, c 0 “ 2 (r 2 ) 1 ^ 2 
have been determined from the equilibrium spectral values. 

In the gravity wave range the integral function in this last com- 
ponent becomes, after using (60) and (65) 

es2 n n ji’ 3dn i ’ (86) 

Thus, for gravity waves after integrating and noting (1) and (66), 
V*(P) * ± (c/c 2 )(e/2) ^ (g/flj) 


’ 1 n 1 k i X(k i* n i )dk i dn i 1 (c/c o } k 
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■ ± (6/2) 1 / 2 (n/n l ) . ( 87 ) 

Combining (87) with (67), 

c 2 ■■ c 2 [l + B(n/ ni ) ± (a/ 2 )i/ 2 ( n / ni )] . ( 88 ) 


Since 6 is such a small quantity, the square root of e is an order of 
magnitude greater than 6. It is seen from (88), then, that the random 
scattering in the gravity wave range has a far greater effect on the 
phase speed than the mean deviation. Physically, this large scattering 
arises from quasi -Euleri an velocity components which have a small mean 
but large fluctuations. 

In terms of wave number and frequency, (1) and (67) are substituted 
into (88) to yield 

n 2 . gk[l + e(n/n 1 ) ± (e/2) 1 / 2 (n/n 1 )] . . (89) 

Once again using a truncated binomial series, 

gk = n 2 [l - 6(n/n 1 > ± ( B/2) 1 / 2 (n/n l ) ] . (90) 

As in the expression for the phase speed (88), the random scattering has 
a greater effect than the mean deviation. 

In the capillary wave range the integral function in the last com- 
ponent of (80) becomes, after using (62) and (72), 


eV /3 f n n 1 '>/3dn 1 . 
Thus, for capillary waves 


( 91 ) 
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V 2 (W 2 = ± (c/c 2 )Y 1 / 3 (3e"/2) 1 / 2 n 1 1 / 3 , (92) 

Combining (92) with ( 73 ), 

c 2 - c 2 {l - (37c 2 )[(3/2 Y k1n n ; - 3c r l /*n } 1 /*‘] 

± (c/$(3372) 1 /y/ 3 0 1 / 3 ^ . (93) 

It is seen from (93) that the large random fluctuations of the quasi- 
Eulerian velocity components cause the random scattering to have a 
greater effect on the phase speed than the mean deviation for capillary 
waves as well as for gravity waves. 

In terms of wave number and frequency, (1) and (72) are substituted 
into (93) to yield, after rearranging, 

yk 3 [l - (3/2)fi'ln + Y ( 1 / 3 )k{[3e' ± (36'/2) 1 / 2 ]nn 1 ( 1 / 3 )} 

- n 2 = 0 . ( 94 ) 

Equation (94) differs from (77) only in the coefficient of y 1 / 3 k. Fol- 
lowing the development of (83), it is readily seen that the cube of the 
solution of (94) for y^k is 

yk 3 = n 2 {l + (3$72)ln n 2 - (n^n) 1 / 3 [3b" ± (3&72) 1 / 2 ]} . (95) 

As in the expression for the phase speed (93), the random scattering has 
a greater effect than the mean deviation. 

The non-linear effects on the phase speed for the gravity and capil- 
lary wave range is shown in Figure 6 and Figure 7. The ratio c 2 /c 2 is 
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random scattering ( — — ) 



A3 


is plotted against n/nj for gravity waves in Figure 6 and for capillary 
waves in Figure 7. From (88), in the gravity wave range 

Cq/c 2 - [1 + gfn/rij) ± (g/2)i/2(n/n 1 )J-l 

* 1 - Bfn/nj) ± (e/2) 1 / 2 (n/n 1 ) (96) 

after using the truncated binomial series. In the capillary wave range, 
multiplying (76) by k 2 gives 

Yk 3 = c^ 2 . (97) 

Applying (97) and (1) to (95) and rearranging, 

e 2 /c 2 = 1 + (3s'/2)1n n : - (rtj/nJVsCSe' ± (H'jZ) 1 /*] . (98) 

In terms of the relationship between wave number and frequency, 
Figure 8 is a gravity wave representation and Figure 9 is a capillary 
wave representation of n 2 /n 2 against n/n : where n 2 is the linear squared 
frequency. For gravity waves 

% = 9k (99) 

and for capillary waves 

n 0 * Yk 3 . (100) 

Figure 8 is determined from (90) and Figure 9 is determined from (95). 

The mean deviation as shown in Figure 6 and Figure 8 decreases in 
the gravity wave range as n/n : increases. In the capillary wave range 
the mean deviation as shown in Figure 7 and Figure 9 increases as n/n 1 
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n n 

Figure 8. Plot of (> /n 2 vs. /n j for gravity waves showing the mean ( } 

and random scattering ( ) 
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Increases* This reversal can be explained by looking at the integrated 
mean quasi -Euleri an velocity in (59)* This term is directly proportional 
to n and inversely proportional to n x in the gravity wave range as shown 
in (67). For capillary waves, however, this term is directly proportional 
to both n and n 1 1 / 3 as shown in (77). A change in n/n^^ will therefore 
cause the velocity in the two wave ranges to change in opposite direc- 
tions. 


The random scattering for gravity waves as depicted in Figure 6 
and Figure 8 increases as n/n l increases. For capillary waves, however, 
depicted in Figure 7 and Figure 9, the scatter decreases as n/n x increases. 
These results are explained by considering the integrated quasi -Euleri an 
velocity spectra from which the random errors are derived. In the grav- 


ity wave range the spectrum is computed from the square root of the in- 


tegral function in (58) after using (60). Thus, 


"^(k^n^dk^] 1 / 2 = (e/2) 1 / 2 (g/n j ) . (101) 

Equation (lo'l) shows that the velocity increases as decreases. For 
a constant n, therefore, the velocity will increase as the ratio n/n x 
Increases. This causes the increasing scatter in Figure 6 and Figure 8. 
For capillary waves, (58) and (62) yield 


[ k n "i x <ki» n 1 >*i d n 1 ] 1/2 = y (102) 

so that the velocity decreases as n 2 decreases. Thus, for a constant 
n, the scatter will decrease as n/n 1 increases as reflected in Figure 7 
and Figure 9. 
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Appendix A 

Third Order Relationship Between dA(k,n) and dB(k,n) 
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Rewriting (31) 

-indB(k,n)exp[i(k-x - nt)] = |k|dA(k,n)exp[|kJ c]exp[i(k-x - nt)] 

+ ^ n k-k i dA(k ,n)dB( k x ,n 1 )exp£ | k | ^]exp{ i [(k+k 1 ) -x - (n+n^t]} 

(31) 

Substituting a series expansion for exp[|k|?], and using (26), 
j^-indB(k,n)exp[i (k*x-nt)] = j^|k| dA(k,n)[l + |k|c + j|k| 2 ? 2 ]exp[i(k*x-nt)] 

+ knkjn -*-i dA ^* n ^ dB ^i» n i^ 1 + Ms! c3 e xp 

{iKk+k^’X - (n+n^t]} 

= ^|k|dA(k,n)exp[i(k*x - nt)] + (nk 

exp{i[(k+k 1 ) *x - (n+n^t] + \ (nk n k n ! !S ! 3 dA( k ,n)dB ( k 2 

- -l 1*2 2 

dB(k 2 ,n 2 )exp{i[(k+k 1 +k 2 )-x - (n+n 2 +n 2 )t] + ^ ^ k*^ 
dACk.nJdBfk^n^expO’Kk+k^*^: - (n+n^t]} 

+ k^ 1 ^yj^!5-!5i dA( ^ n ) dB ^i» n ^ clB ( k 2 » n 2 ) 


expdEfk+kj+^J-x - (n+n 1 +n 2 )t]} 


(A-l) 
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For the first order let 

k = k Q , n = n Q . (A-2) 

For the second order let 

k + kj = k Q * n + n i = n 0 * { A-3) 

For the third order let 

- + ~i + -2 = ~o ’ n + n i + n 2 = V * (A 1 - 4) 

Then, after substituting (A-2), (A-3) and (A-4) into (A-l), 


/ / 


k n - in o dB ^ 0 * n o )exp[i ^o'?- n o t)] * k n dA( VV exp[i(l V*- n o t)] 

.n -Ci n 


0 0 


+ k n ic n I V^il 2dA( -o'-i* n o" n i^ dB ^i* n i^ exp f i ^o'?" n o t ^ 
-0 0-1 1 


♦iklkiy, l -o*-r-2 13dA{ -o'^r l 52»v n r n 2 ) 

-0 0-1 1-2 2 


dB(k 1 ,n 1 )dB(k 2 ,n 2 )exp[i(k 0 -x-n 0 t)] + n^rij 


dA (V^*V n i )dB( ^‘ n i )exp[i(!5 o^- n o t)] + 


l^r^l ( ^r^ ),, 5i dA{ V^-^ > v n r n 2 ) 


(A-5) 


dB(k x *n x )dB(k 2 ,n 2 )exp[i(k 0 *x-n 0 t) ] . 
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Since each component of (A-5) Is Integrated over ( k Q »n Q ) » only the Inte- 
grands of the Integration over ( » n 0 ) need to be considered. After 
dropping the zero subscripts, these integrands in (A-5) can be rearranged 
so that 


dA(k.n) = - dB(k,n) - ' ' [IHJ 2 + !' 1 *(H,)3 "(k-kj.n-n, 


yjjy tlA (^‘!S 1 *!$ 2 * n ' n 1 ' n z ) c|B (!< 1 »n I )dB(k 2 ,n a ) 


(A-6) 


Using successive approximations, to the first order 


dA(k,n) « - yj~j- dB(k,n) 


(A-7) 


To the second order, after substituting (A-7) for dA{k-k i .n-n^ in (A-6), 
dA(k,n) - - ^ydB(k,n) - ^[Ik-kJ 2 + [. 

dBfk-kj.n-rijJJdBfkj .n^ 


In 


. . ./ /' MHi) 

M + K „ tt i - ivJ 


■jYp dB(k,n) + fg -ffk'-k I (n-njJdBCk-kj.n-nj) 


dBfkj.n^ 


(A-8) 



Continuing the approximation, substitution, of (A-8) into (A-6) yields 
the third order relationship 


dA(k.n) = - 


in r r k* ( k-k ) 

-^ r dB ( k , n ) . y 


Ifn - nj ) 

jpg- r dB(k-k 1 ,n-n 1 ) 


. / 1 

+ 1 k 2 n 2 TF k l l l k ‘ k 1 ‘ k 2 l (n - n r n 2 )dB (Hr!s 2 * 


n-n 1 -n 2 )dB(k 2 ,n 2 ) 




+ ^~~1~~2 ^ 1 _ 1^ JFf 


L 


i(n-n 2 -n 2 ) 
J k-k -k I 

1 - >i -i 1 


dB ^^r!s2* n -‘ n r n 2 ) 


dB(k 1 ,n 1 )dB(k ,n ) . (A-9) 


Rearranging 


dA(k,n) = - jp- dB{k,n) + i k 1 nJk][k-kJ. (n ' n i )dB( H 1 * n - n 1 ) 

r r r r (n-n -n_) 
dB <‘iV + 1 k.n.k.n. [kj ~ 


- 1 . 1-2 2 

k-Ck-kiJCk-kjHk-kj-k.,) 

I k _ k ll l k - k !- k 2 l ’ 


j ( k + k r k 2 )-{ k - k r k 2 ) 


dB ( k ' k i - k 2 * n ' n r n 2 ) 


dB(k 1 ,ri 1 )dB(k 2 ,n 2 ) 


(A-10) 



Appendix B 

Third Order Expansion of the Dispersive Relationship 
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Rewriting (39) 

95 = - {3*/3t) c - + Y (vgc)[1 + (v h? ) 2 ]' (3/ - 2) . (39) 

Expanding the term raised to the -3/2 power in a binomial series, 

95 ■ - ( 3<J>/ 3t> c - ^(v*) 2 + Y (v 2 c)[l - (3/2)(v h ?) 2 ] . (B-1) 

Using (25) and (26) to substitute for $ and 

kn 9dB(k,n)exp[i(k'X - nt)] = indA(k,n)exp[|k|i;]exp[i(k-x - nt)] 

" S' kn^n^ M sM!5i 1 - k-k 1 )dA(k.n)dA(k 1 .n l )exp[(|k| + |kj)c] 

exp{i[(k+k,).x - (n+n^t]} - Y " jk| 2 dB(k,n)exp[i(k-x - nt)] 

- Y(3/2) Hiiii |k| 2 k 1 -k 2 dB(k > n)dB(k 1 .n.)dB(k,,n,) 

- -1 1-2 2 

expdEtk+k^J-x - (n+n +n )t]} . (B-2) 

Expanding the surface slope exponential functions into Taylor series 
around s=0, 

^gdB(k,n)exp[i(k.x-nt)] = £' in[l+|k| c + j|kj 2 c 2 ]dA(k,n) 
exp[i(k-x - ntJj-l^^lkllkJ-k.kjJCHIkl+lkJJc] 



58 


dAO^nJdAfkj ,n 1 )exp{i[(k+k 1 )«x (n+n^t]} 


- |k| 2 dB(k,n)exp[i(k-x-nt)] - y(3/2) l^k n k n 


- -1 1-2 2 


dBfkjnJdBfkj »n 1 )dB( k 2 ,n 2 )exp{i[(k+k 1 +k 2 ) -x - (n+n x +n 2 )t]} . 


(B-3) 


Rearranging, 


£n 9^B( k ,n)exp[i ( k -x-nt) ] - ^ indA(k,n)exp[i(k*x-nt)] + £££ ^ 

1n|k|dA(k,n)dB(k i ,n 1 )exp{i[(k+k i )-x - (n+n 1 )t]}+ 1 ^ £ £ £ 

in|k[ 2 dA(k,n)dB(k 1 ,n 1 )dB(k 2S n 2 )exp{i[(k+k 1 +k 2 )*x - (n+n 1 +n 2 )t]} 

/ 

“ F knk n ~ k*!^) ^A(k J n)dA(k i ,n 1 )exp{i[(k+k 1 )*x - fn+n^t]} 

• J knk^^^l-l l-il " + I !$ ! I ) dA(k ,n) dA( k 2 ,n j ) 

rfB(k 2 ,n 2 )exp{i[(k+k 1 +k 2 )-x - (n+n 1 +n 2 )t]> - y ^ |k| 2 
dB(k,n)exp[i(k*x-nt)] - y(3/2) f n [ J n 1^1% -k 2 dB(k,n) 
dB(k 1 ,n 1 )dB{k 2 ,n 2 )exp{i[(k+k l +k 2 )-x - (n+n^n^t]} . (B-4) 
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Numbering each term on the right hand side of (B-4) consecutively from 
I to VII and substituting for dA(k,n) from (36), to third order 


I = 


kn ]fr dB (k,n)exp[i{k-x - nt)] - n(n - n t ) 

rrr r r r n(n-P -n„) 

dBCk-kj ,n-n :L )dB(k i .n^expCifk-x - nt)J - ' * * * * ' 1 


knk 1 n 1 k 2 n 2 


TFT 


1 k*(k-k, )(k-k. ) * 

7<HrV*IHr!5 2 l- n' Hi l lH x -k 


k-(k-k 1 )(k-k 1 )-(k-k 1 -k 2 ) 

rr^"~ ’ 


dB ^'^-!52* n - n r n 2 )dB ^i * n 1 )^B(k 2> n 2 )exp[i(k.x -nt)] . (B-5) 


For second order let 


k - ki = r, n - n 2 ?=.. n" . (B-6) 

For third order let 

k “ kj - «V. n - n i - n 2 = n" . (B-7) 

Then, after substituting (B-6) and (B-7) into (B-5) and dropping the 
primes, 

r r n 2 rrr . k*(k+k) 

1 = kn TIT dB (k.n)exp[i(k-x - nt)] - knk^ J ' k | fk^ | "(" + V 

dB(k,n)dB(k 1 ,n 1 )exp{i[(k+k 1 ) -x - (n+n^t]} 

////// (n+n i+ n 2 > n r 1 , , (>S + k 1 + fS 2 )-(k+k 2 ) 

- n -i n i-2 n 2^ + -i +k 2l L 2 * - 1 ~ " |k+k 2 |jk| 
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(k+k a ) «k]dB(k v n)dB(k 1 »n 1 )dB(k 2 ,n 2 )exp{it(k+k 1 +k 2 ) 


- (n+n 1 +n 2 )t]> 


(B-8) 

Looking at the second term on the right hand side of (B-4), after using 
(36) 


II * n 2 dB(k,n)dB(k i> n 1 )exp{i[(k+k i )-x - (n+n^t]} 


- -l l 


iff f r r 

knk i n i k 2 n 2 lF^T n( """i )dB( -'^i » n i> 


dB(k 2 ,n 2 )exp{i[(k+k 2 )-x - . (n+n 2 )t]} . (B-9) 

For the third order let 

k - kj = k' , n - rij = n" . (B-10) 

Then, after substituting (B-10) into (B-9) and dropping the primes. 


11 = knk n n 2 dB(k,n)dB(k 1 ,n l )exp{i[(k+k 1 )-x - (n+n^t]} 
rff f t r (!$"^bi)‘!s 

" tok h k i — TB1 n(n + n 1 )dB(k,n)dB(k l> n 1 )dB(k 2 ,n 2 ) 


exp{i E(k+k 1 +k 2 ) -x - (n+n J +n 2 )t]} 


(B-11) 



Using (36), the third terra on the right hand side of (B-4) becomes 



fit S S S 2 
knk 1 n 1 k 2 n 2 


fkJdB(k,n)dB(k i 


,n )dB{k ,n )exp{i[(k+k ,+k )• 

1 ~Z 2 - ~ 1 '2 


- (n+nj+n^t;)} 


and the fourth term becomes 


(B-12) 


IV = - 1 / 


2 knk^HkJIkl- k.kjj - dB(k.n) + 1 ' ' 

v 2 2 ** -2 


(n-n 2 )dB(k-k 2 ,n-n 2 )dB(k 2 ,n 2 )| J - 3- dB(k,. ni ) 


expdCfk+k^-x - (n+n )t]} . 
Rearranging, 


(B- 13) 


IV = 1 fff f 
2 knk 1 n i 


1 - 


Hi 

HTTFT 


nn^Bfk.nJdBCkj.njJexpHCfk+k^.x 


- (n+n^t]} - 


1 111 I 1 S 

2 knk J n l k 2 n 2 


i - r '- 1 
. 7 ¥TT ¥7 




n(n 1 -n 2 )dB(k,n)dB(k 1 -k 2 ,n i; n 2 )dB(k 2 ,n 2 )exp{i['(k+k ).; 


-(n+n )t] - i S!S 1 l 1 
1 2 knk n k n 

- 1 1 o 


■ m r\ II 

~ -1 1-2 2 


k-k 2 

TITltj 


.MHO 

TFy~ 


n /n-n 2 ) 


1 
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dB(k-k 2> n-n 2 )dB(k 1 ,n 1 )dB(k 2> n 2 )exp{i[{k+k i )*x - (n+n^t]} 


(B-14) 


In the third term on the right side of (B-14), let 


k - k 2 = k", n - n 2 = n" . 


(B-15) 


Then, after substituting (B-15) into (B-14) and dropping the primes, 

k«k 


IV = 1 M * 

2 knk 1 n 1 


1 - 


- -i 


"TTT17" 


nn 1 dB(k,n)dB(k 1 ,n z ) exp{i [(k+k^ *x 


- <*«!>« ~ 7 knk n k n 1 ' 
- -1 1-2 2 


k-k 






n(n i -n 2 )dB(k,n)dB(k 1 -k 2 ,n 1 -n 2 )dB(k 2 -,n 2 )exp{i[(k+k i ) *x 


(n+n x )t]) ' 1 kik^k^ 


■ fk+y*! 


(k+k 2 )-fe 

' FT 


(nn l )dB(k,n)dB(k l ,n 1 )dB{k 2 ,n 2 )exp{i[(k+k l +k 2 )*x 


- (n+n 1 +n 2 )t]} . 


■(B-16) 


Substituting (36) into the fifth term on the right hand side of 
(B-4) and rearranging. 


7 taM.M ClkllU * HJdkl+lkJ) 


-- 11-22 


v-m 


ln dB(k.n) 



dB(k ,n )j 


dB (k 2 » n 2 )expfi E(k+k a +k 2 ) *x 


(n+n,+njt]) 
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e l ////// 

2 l<nk x n x k 2 n a 1 - 

dB(k 2 ,n 2 )exp{i[(k+k 1 +k 2 )-x - (n+n 1 +n 2 )t]} . (B-17) 

The dynamic boundary condition is now specified in terms of dB(k»n) so 
that 

gdB(k »n)exp[i (k-x-nt)] ® I + II + III + IV + V + VI + VII * 

(B-18) 

Equation (B-18) can be rearranged and then rewritten as 

F(k,n; k 1 ,n l ; k 2 ,n,,)dB(k,n)exp[i (k-x - nt)] = 0 (B-19) 


(|k| +|k 1 |)nn 1 dB(k,n)dB(k l ,n 1 ) 



ro|— » 



64 


d B (k 1 ,n l )dB(k 2 ,n 2 )exp{i[(k 1 +k 2 )-x - ( 0 j+n 2 )t]} 


1 / / / / 

F k,n,k„n„ 
-1 1-2 2 


1 - 


k - k i 


k [| k 1 ” 


k j • ( k x -k ) 

j lc p k " 2 | " (n i~ n 2 )dB( -l*~2* n i' n 2 > 


dB(k 2 ,n 2 )exp[i(ki*x - r^t)] 
Since dB(k,n) can be any function. 


(B-20) 


F(k,n; k 2 ,n 2 ) = 0 


(B-21) 


Thus, after setting (B-20) to zero and rearranging terms, the third 
order dispersive relationship can be written in the following form: 

, n 2 rr Mk+kj) 1 k-k. = 

9 + Y l-I 2 ‘ w = M, n ' TFTfFrr n(n+n i ) + 2 (1 ■ iFHkT )nn i 

f , t r fn{n+n +n ) 

dB(k 1 ,n 1 )exp[i(k 1 -x - n^)] - 

1 /ll01 x , (!? + !5i + k 2 )-(k + k 2 )(k + k 2 )-k 

F (*5 +2 !< 1 ) *!< [k+iTTlkl 


- -2 


(k+k )«k , 1 

+ | k t ’ n{n+n i) - F n2 ^l - F 


1 - 


TFTT¥7 


(|k|+|k 1 |)nn 1 


+ F 


1 - 


(k+k 2 )'k 


‘WTrnr 

- -2 - 1 


]T| nn ! + 2 v|k| 2 k^kj 


dB(k 1 ,n 1 )dB(k 2 ,n 2 )exp{i[(k l +k 2 )-x - (n 1 +n 2 )t]} 
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ISIS 



( *ry 



n(n 1 -n 2 )dB{k 1 -k 2 , lyn^ 


dB(k 2 ,n 2 )exp[i(k 1 -x - rijt)] 


(B-22) 
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2* Statistical properties of surface waves under non-uniform currents. 

2.1 Introduction 

Previous studies (Longuet-Higgins and Stewart, 1961; Huang et al, 

1972; Phillips,. 1966) have shown that when waves encounter current, 
wave characteristics undergo changes due to interactions between the 
waves and the current. 

In order to utilize the phenomenon of wave-current Interactions as a 
means of measuring current, it is desirable to first conduct a comprehensive 
study of the effects of current on relevant statistical properties of 
waves. This includes those of wave elevation, wave amplitude (peak), 
zero crossing rate, number of maxima, velocity of zeros and specular 
points, and other related quantities for both the one-dimensional and two- 
dimensional wave systems. Those quantities that are most sensitive to 
wave-current interactions and particularly suited for current measurement 
by remote sensing devices will then be identified. 

Work to date has been devoted to one dimensional wave systems only. 

The effects of current on wave elevation, wave amplitude, zero crossing 
rate, number of wave maxima, and velocity of zeros have been studied and 
reported in the following. • 


2.2 Distribution of Wave Elevation 
Both theoretical and field studies show that the surface elevation of 
a random wave field is approximately Gaussian. That is, the probability 
density function p(0 of surface elevation £, measured from mean sea 
level is 


V5na 


exp[- 


ii!] 

2 „2 J 


p(0 = 


( 2 . 1 ) 
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In which cr is the standard deviation of £ which can be obtained from the 
spectrum of surface elevation. 

It was shown previously (Huang et al, 1972) that under the influence 
of a steady non-uniform current, the frequency spectrum of a random 
gravity wave field is given by 


<f>(n) = 


[1 + (1 + ^) 1/2 ] [(1 + -^ 2-) 1/2 + (l + ^)] 

8 6 S 


in which n is total frequency, U is current speed, g is gravitational 
acceleration and <J>*(n) is wave frequency spectrum without the influence 
of current. In this study, the Kitaigorodskii-Pierson-Moskowitz spectrum 
is used. That is, 

4>*(n) = exp[-a(% 4 ] (2.3) 

n * 5 

_ o 

in which a and 3 are non-dimensional constants equal to 0.74 and 8,1 x 10 
respectively, and n Q = g/W, W being mean wind speed. 

The standard deviation a that appears in equation (2.1) is given by 
' n c 

a ® f 4>(ti)dn (2.4) 

o 

in which n c is the cut-off' frequency of the wave spectrum taken as that 
of a wave 30 cm in length. 

The probability density functions as given by equation (2.1) are 
plotted in figures 2.1, 2.2, and 2.3 for different current and wind 
speeds. It is seen that for a given wind, positive current reduces the 
probability of large wave elevations while negative current increases the 
same. It is also noted that the larger the wind speed, the smaller is 
the effect of current on wave elevation. 



U ** 4 tn/sec. 



Figure 2.1 Probability density fu 
under current for wind 






1-D current 



Figure 2.2 Probability density function of wave surface elevation 
under current for wind speed at 10 m/sec. 


1-D current 



Figure 2.3 Probability density function of wave surface elevation 
under current for wind speed at 20 m/sec. 
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2.3 Distribution of Wave Amplitude 
For a Gaussian, zero mean, stationary sea, the probability density 
function p(n) of the amplitudes rj (peak) of the waves is given by 
(Longuet-Higgins , 1956). 

2 


P(n) - 


V2 


inn. 


[ e exp (- \ ^ 2 ) + 

* T&O £ 




^7 n exp( - * ill ' 



exp(- < 2 * x 2 )dx] 


(2.5) 


in which 


e 2 = 



m o m 4 


( 2 . 6 ) 


a parameter measuring the r.m.s. width of the spectrum and, in the case 
when the waves are considered a function of time t, 

m 0 = J 4>(n)dn 
n 

m 2 = X n 2 <f>(n) dn (2.7) 

n 

m 4 ” X ( n ) dn . 

n 


If the waves are considered as a function of space n, 4>(n) is to be 
replaced by ^(k), the wave-number spectrum and n is replaced by k, the 
wave number. In either case the integrations in equation (2.7) are to 
be extended over the entire gravity wave range. 

Other statistical properties of wave amplitudes such as the mean 

» __ 

Uj_, variance i> 2 » and coefficient of skewness B, can be readily derived 

from the expression of probability density function p(n) of wave peaks 

in equation (2 .5) . These are 





72 

U 1 = >J-5 m o(l- e2 ) , 
* 


(2.8) 

U 2 88 in 0 [1- (|* -1) (1- e 2 ) 3 » 


(2.9) 



(2.10) 


Another quantity of interest is the proportion of negative peaks (peaks 
of negative magnitude) out of the total number of peaks of the waves, 
denoted by y, and is given by 

Y - j [l-(l-e?) 1/2 ] (2.11) 

Effects of current on the statistical quantities mentioned above are 
computed for the case when the waves are treated as functions of time and 

presented in graphical forms. In figure 2.4 the ratio of mean wave 

» 

amplitude with and without the influence of current is plotted against 
current speed U with mean wind speed W as a parameter. It is seen that 
positive current reduces mean wave amplitude while, due to energy pile up, 
negative current increases the same and the effects of current are more 
pronounced at lower mean wind speed. 

Figure 2.5 gives the ratio of variance of wave amplitudes u 2 with 
and without current. The same trend that is observed in figure 2.4 is 
noted here. 

Presented in figure 2.6 is the ratio of coefficient of skewness B 
with and without current considered. It is seen that positive (negative) 
current gives rise to an increase (decrease) in skewness of wave amplitude 
distribution. This is because under the influence of a positive current 
the waves become smoother and therefore of narrower band giving rise to 
more positive peaks than negative peaks resulting in a more skewed 





-0.5 

4 

0 0.5 1.0 

Figure 2.5 Change of variance of wave amplitude under current. 
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amplitude distribution. Conversely, negative current feeds energy into 
the wave system. The waves become more choppy and therefore of wider band 
giving rise to a less skewed wave amplitude distribution. 

Figure 2.7 presents the ratio of y with and without current included. 
That there are less negative peaks in the presence of a positive current 
than when the current speed is negative is clearly noted. The physical 
mechanism underlying the behavior of the curves in figure 2.6 also governs 
that of the curves shown in figure 2.7. 

With the statistical moments of wave amplitudes discussed, the 
characteristics of the probability density function of wave amplitudes as 
given by equation (2.5) and presented in figure 2.8 can be readily under- 
stood. In figure 2.8 the probability density functions p(n) of wave 
amplitudes are given for wind speed W=10 m/s and current speed U=2, 0, 
and -0.4 m/s. That positive current gives rise to smaller values of mean 

and spread of wave amplitudes than negative current does is clearly seen. 

... 

The curves also show that the amplitude distribution is more skewed under 
positive current than in the case of negative current. 


2.4 Expected Number of Zero Crossings and Maxima 
In addition to amplitude distribution, the statistics of number of 
threshold crossings and extrema of random functions and, in this case, 
sea waves are often used. While it is difficult to obtain the distri- 
butions of these quantities, their expected values are easy to compute 
both theoretically and from field observations. 

According to Longuet-Higgins (1962) , the expected rate of zero 
crossings, denoted N Q , is given by 

N = - A 1 ' 2 

o ir x m n 


( 2 . 12 ) 
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and the expected number of extrema per unit time is 

H e - 7 ' < 2 - 13 > 

When the waves are considered as a function of time, the m f s are obtained 
from equation (2.7). When the waves are treated as a function of space, 
4>(k) and k should be used in place of <fi(n)dn in equation (2.7). In this 
study, only the former case is considered. . 

In figure 2.9 the ratio of N 0 with .and without the influence of 
current is presented. In the presence of positive (negative) current, 
waves are of a narrower (broader) band resulting in a decrease (increase) 
in expected number of extrema, per unit time. 

In figure 2.10 the ratio of N 0 with and without the influence of 
current is also presented. It is seen that the same trend observed in 
figure 2.9 is noted. 


2.5 Velocity of Zero 

The quantities examined above all refer to waves either treated as 
a function of time t for a fixed point in space or as a function of space 
x, for a given instant of time. Moving waves are, however, functions of 
both space and time. It is, therefore, of interest to investigate some 
statistical properties of moving wave forms. Presented in the following 
is velocity of zero. 

Consider a one-dimensional random wave f(x,t), the velocity of the 
zero, denoted C, is given by (Longuet-Higgins , 1956, 1957) 


3f/5t 
3f /3x 


(2.14) 


The distribution of C, it was shown (Longuet-Higgins, 1956, 1957), depends 
on t|#(k), the wave-number spectrum. Under the. influence of current. 



Ne/(Ne) o 

- It 3 

- 1.2 

- 1.1 


1 



81 


ij<(k) is given by (Huang, et al, 1972) 

♦ 

= ^^7 1? exp[ -“ ( k(fe) )4] <2 - 15) 

in which c = (g/k) i 2 / 2 , the phase speed. The above expression for $(k) 
is derived based on the assumption that when there is no current the wave 
frequency spectrum is given by equation (2.2). Defining the following 
moments of ij>(k) as 


m 2 = J* k^t|>(k)dk 
k 

f kn^(k)dk (2.16) 

k 

m o ~ / n 2 4>(k)dk 
k 

in which under the influence of current the relationship between frequency 
n and wave-number k is given by (Huang et al, 1972) 


gk 


n‘ 


[! + i ( 1 + ir )1/2]2 


(2.17) 


the probability density function p(£) of velocity of the zero C is 
(Longuet-Higgins , 1957) 


p(C) 


V m 2 


[(c-o 2 + A 0 /m 2 ] 3/2 


(2.18) 


i » 

in which A 0 = m 2 m'^ - mj, z and £ - -m^/u^, the mean of £. The interquartile 
range td, a measure of the spread of 5, is m = 2%/3A 0 /m2. 

Plotted in figure 2.11 is the ratio of If with and without current. 

It is seen that the mean of £ increases (decreases) with increasing positive 
(negative) current. It is also noted that the larger the wind speed, the 
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less Is the effect of current on the quantity considered. It is worth 
noting that comparison of figure 2.11 with figures 2.4 to 2.10 indicates 
that the effect of current is more pronounced on velocity of the zero 
than on any of the other quantities examined in this report, suggesting 
that velocity of zero may well be a quantity to be used for current 
measurement by remote sensing devices. 

In figure 2.12 the interquartile range to with and without current is 
studied. The curves exhibit the same characteristics as those of the 
curves in figure 2.11. 

Finally, the probability density function p(?) of £ is shown in 
figure 2.13 for W = 10 m/s and U = 2 m/s, 0, and -0.4 m/s. The influence 
of current on the mean and spread of velocity of zero shown in figures 
2.11 and 2.12 is clearly reflected in figure 2.13. 



-0.5 


0.5 


1.0 


1.5 


2.0 


Figure 2.11 Change of phase velocity at zero crossing points . 
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Figure 2.13 Probability density 
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Conclusion 

The results of the theoretical studies indicate that the statistical 
properties of the sea surface are sensitive to whatever causes changes 
in the wave energy distribution; i.e., the spectral functions. This 
is true even with higher order expensions such as in the Gram-Char lier 
series for wave height probability density function. In the open ocean, 
mechanisms that will change energy distribution substantially are, of 
course, winds and major currents. However, before more detailed knowledge 
of generation of waves by wind is available, It is impossible to utilize 
these properties as a means to measure surface wind over the ocean. On 
the other hand, current-wave interaction is more definite and the effect 
is of the first order. Furthermore, since the change depends on relative 
values of currents , it is less susceptible to the error in establishing 
an absolute relationship as required in wind wave generation studies. 

Under normal conditions a wind system will cover an area substantially 
larger than the wave scales and thus provide a homogeneous field for wave 
actions. Take an inlet or a river mouth for example, the local flow will 
generate a non-uniform velocity field over a homogeneous wave field. 
Consequently, the waves will interact and change characteristics depending 
on the flow conditions. Such changes can be detected easily from the 
probability density function. 

Another example is the Gulf Stream. The width of the system is about 
100km with strong shear zones on its sides. Though the decay scale of a 
wave 30cm in length is only 1700m, it is not unusual to find substantial 
current velocity changes across such a distance as reported by Stommel 
(1966). When waves encounter such current systems, their characteristics 
will change accordingly. Although the analyses in this report did not 
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include shear current cases where as in the Gulf Stream the shear will 
be important, a qualitative prediction on the wave height probability 
density function can be made with the available results as a guide. To 
facilitate the discussion, a schematic diagram is given in figure 4.1 in 
which a typical velocity cross section of the Gulf Stream measured by 
Worthington (1954) is shown. The total flow region is further divided 
into four different sub-zones numbered one through four, with A and B 
indicating the local extreme of velocity. If wind generated waves are 
propagating from the open ocean of zone 4 across the Gulf Stream to 
zone 1 toward the land, they will first encounter counter current in 
zone 3. The probability density function of wave elevation will become 
increasingly flat until it reaches point B. From there on the relative 
change of the current is increasing in the direction of the waves; there- 
fore, the flatness of the density function will become less and less until 
point A is reached. After waves pass point A, relative change of the 
current is in the opposing direction of the wave again, and the flatness 
will increase accordingly. Granted that no single wave will last the 
whole current system under moderate or light wind conditions, but as long 
as the locally generated waves can experience different current conditions, 
the changes in probability curves will still show the effects. Such 
changes are easily detected by means of remote sensing devices, and thus 
it offers a new way of measuring the major ocean current systems. 

Although this seems to solve the problem in principle, there are 
related problems that will have to be solved before the final result can 
be used to set up a routine in practice. One of the problems is that the 
results in the theoretical analysis are based on one— dimensional currents 
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Figure 2.14 Schematic diagram on change of wave elevation distribution 
the Gulf Stream. 


across 


89 


and waves, which is not realistic compared to all the ocean situation; 
therefore, an extension into the two-dimensional model is necessary* 

Secondly, we have to transform all the analyses into wave number 
space. For the one-dimensional model, a frequency spectrum is sufficient, 
while frequency spectra are easily available from data collected by the 
conventional methods in the form of time series, the remote sensing method 
is more effective in collecting data at any instance covering a large 
area. Such data in spatial variables will give us more complete information 
in the form of wave number spectra; therefore, the results presented in 
this report will have to be converted into parameters reflecting spatial 
distributions. 

Thirdly, we need a definite relationship between wave number and 
frequency in a random wave field. This is critical especially in this 
developing stage when checks between remote sensing and in situ results 
are indispensable in calibration and improvement of the system. In such 
processes, it is necessary to transform information from frequency space 
(obtained by in situ methods) into wave number space (obtained by remote 
sensing methods) and vice versa. In the past a rough relationship based 
on a single train of linear small amplitude waves has been used. The 
deficiency of this approximation is tolerated for lack of alternatives. 

A new method of establishing a more accurate relationship is being 
developed during this study. 
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3. Interaction of Capillary-Gravity Waves with Non-Uniform Current 

One of the problems of a laboratory study of wave interactions has been 
the difficulty of generating waves in the laboratory which correspond with 
the type of waves involved in the existing mathematical descriptions of the 
interaction process. One such description is that of Huang et al (1973) 
which rigorously develops the interaction of pure capillary waves with current 
in the one-dimensional case. However, the single wave trains generated under 
laboratory conditions are almost never pure capillary waves, but are rather 
capillary-gravity waves, described by the equation 

+ gk, (3.1) 

where a is the wave frequency in radians, y is surface tension constant, 

g is the gravitational acceleration constant, and k is the wave number, equal 
to 2^/A, where A is the wavelength. Waves generated with a frequency of 30 
cy/sec still consist of approximately 34% from gravity wave contributions, 
represented by gk in (3.1) above. At 40 cy/sec, these gravity wave contri- 
butions are still approximately 23%; thus increasing the frequency results 
in the generation of waves which approach the pure capillary case. 
Unfortunately though, one phenomena limits the development of a single train 
of pure capillary waves — cross waves generated at the wave maker perpendicular 
to the direction of propagation desired. Between 30 and 40 cy/sec, these 
cross waves appear and increase in intensity until at about 40 cy/sec almost 
all energy goes into these cross waves. At this point, no single wave train 
is possible. Thus in reality, the single wave trains generated in the 
laboratory are basically capillary-gravity waves at the shorter wavelengths, 
and conversely, gravity-capillary waves at the longer wavelengths. Therefore, 
what must be described mathematically is a more complex wave made up of both 
capillary and gravity type terms. In general, phase speed c and wavelength 
A are related as shown in Figure 3.1. To the left of the minimum phase speed, 
at x = 1.76 cm, are the values corresponding to capillary-gravity waves, and 
to the right, gravity-capillary waves. Moving further to the right or left 
results in a more "pure” wave. Note that for any c other' than the minimum 
value, the relation between c and A is double-valued. Because of this, only 
the capillary-gravity region will be treated here. The method to be used is 
the perturbation method, with the perturbation e being such so that e < 1. 




Figure 3.1. The relationship between phase 
speed c and wavelength X 
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In this approach, we use the energy conservation equation as used by 
Phillips (1966) which states that 


3U fi 

H + hr ilLlV * + C S« ]} + S aB taT ° 6 > 

01 


(3.2) 


where 


a, 3=1, 2, corresponding to x^ = x, X 2 = y 


c = group velocity 
“a = current velocity 
E ■ wave energy 

S a p = excess momentum tensor, or Reynolds stress 
- 6 - dissipation. 

For the one— dimensional case under steady state conditions, and no 
dissipation, (3.2) can be written as 


— — {E[U + c ]} + S - 0. 
3x g xx3x 


(3.3) 


Thus expansions in powers of e are needed to substitute into (3.3). We begin 
by rewriting (3.1) as 

♦ o 2 « yh 3 (l + -^7) . 

’ O' 

g 

To develop an expression for c * we take of (3.4) to get 


(3.4) 
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or 


Define now 


2cf~ = 3yk 2 (1 + + YkV 

dK yk - yk 


3yk (1 + - 2g 

8a yk 

— c„ = 


3k * ~ a + -^) 1/2 ' 


(3.5) 


Yk 


K ~ * e ~ ^2 * 

k yx 


where the hat values denote constant values at the points where U - 0. Thus 
using these definitions, (3.5) can he rewritten as 


2 £ o[' 


k 1/2 (i + -y ) 1/2 
K 


3K 3/2 (1 


Ij. * 

. JBa/2 * 


(3.6) 


where c c = \yk , the phase velocity of a pure capillary wave at the location 
where U = 0, thus a constant value. Expanding (3.6) and grouping by orders 
of e results in 


3 - r 

2 C ° L 


k 1/2 -(-4tt) £+ (■ 5 


6K‘ 


3/2' 


24K 


7/2 


) 2 + 




(3.7) 


Now, to develop an expression for U, we start with the kinematic 
conservation equation, as used by Phillips (1966) : 


3k 


8l+ Vn " 


Assuming steady state conditions, this becomes 

7 (a + k*U) = 0* 

For the one— dimensional case along x, this is just 


f^(kU + kc) = 0, 
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since a - kc. This means that 

♦ 

kU + kc = constant. 

A 

The constant will be chosen where U = 0, k = k, c 


c, thus 


kU + kc = kc 


or 


A A 

C C 


1 

K * 


(3.8) 


Since 


= c = 


yk + using (3.1) it follows that 


* yk(l + -£=■) ykK(l + — ■) 

cl = Yk Z = K 

A 2 A p a 

c Yk(l + -fj) Yk(l + e) 
Yk 


or 


K 1/2 (l + -—*) 1/2 
c k 


(1 + e) 


1/2 


Using this result in (3.8) leads to 



K 1/2 (l + 

K 


(1+ e) 


1/2 


A A 1 12 

Expanding, using c = c 0 (l + re - g’E + ••*)* and grouping in terms of e orders 
gives 



K 1/2 )c 


+ !£(!_ 1 — .) E - (1— + I) £ 2 

z° k 3/2 ; 8^° V k 7/2 k' e 



(3.9) 


It thus follows that 




9U 

3x 



2K 


1 % a 

l/2 )c o 


2K 


5/2 


K 


. ls t 

E + 8 c o< 


2K 


9/2 


K 



(3.10) 
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Expanding the energy E now as 


E “ E A + eE n + £ E„ + 
o 1 l 


and the Reynolds stress term as 


S xx ■ TTT? “ (3 + £)( 2 " 2 e + b 2 •" )[E ° + ^1 + + '" ] 


or 


S xx “ f E 0 + ( K - E 0> e + <2 E 2 - E 1 + E o> e + ••• » 
when grouped in orders of e. 

By substituting (3*7), (3.9), (3.10) , (3.11) , and (3.12) into (3.3) : 
carrying out the multiplications and grouping by orders of e, the energy 
equation (3.3) becomes, to order t° , 


L| E o[ ( K - 


K 1/2 )c q + |r 1/2 c 0 | [ + |e o c 0 ( _ ^2 


-I! 


1 -1-^=0 


K 2K 


1 / 2 ' Sx 


_± , i 


J. COUJLL.O 


.Multiplying (3.13) by ^. 3/2 
quotient, so that (3.13) can be written as 


*. ^ ^ ^ n A A — 1 J *n.f n 

l.o i.n a pai L4.ai. -l. v <x u j- *■ ^ x ^ 


J_ E 0 Cq[k + ^ ] 
3x „3/2 


= 0. 


K 


Thus 


E o £ o (| + |k 1/2 ] 

„3/2 


= constant. 


Again we choose the constant at the locations where U — 0, for there 
E 0 = E 0 , K = K = 1. This results in a. constant of 


E o c o ^ + 2 ^ 3£ a 

= constant = 2 E o c o * 


(3.11) 


(3.12) 


(3.13) 


(3.14) 


(3.15) 


Equation (3.15) now becomes 



or 


*o°o'k + I r1/2] 3. . 

3/2 ^o c o » 

* K 


E, 


3^3/2 

7T 


E 0 1 + V/2 • 

k r 


*6 


or 


(3.140 


Equation (3.16), where the order is e°, corresponds to the pure capillary 
wave case. This result is the same as that obtained by Crapper (1972), 

Huang et al (1973), and Holliday (1973). Using (3.9) and (3.10a), (3.16) may 
be rewritten in the form 


4-(£-) 3 


U 


3 ,c 


(3.17) 


r + 
o 


Cfi 2 Cq 


Writing the energy equation (3.3) in its expanded form to order e , we 
obtain 


3 „ * 7 l 

3x E o c o^2K 


4K 


3 , , _ a 9 , 1 , 1 > 3K 

3/2 } + E ° C ° ( - OT ,5/2 + /tr 2 + nT ,l/2 } 3x 


8K' 


4K 2K 


.3 „ a ,1 . 1„1/2 X „ a , 3 , 3 v 3K _ „ 

+ t — E-, c (— + ipi ) — E. c ( — ;r + — TTo) “ 0 • 
3x 1 o K 2 ' 1 o^„2 ,„l/2' 3x 


(3.18) 


2K AK J 

Note that the E-i terms of (3.18) are of the same form as the order e° case. 
Thus (3.18) can be simplified immediately by multiplying through with 


K 3/2 ’ 


resulting in 


1 !_ E ; , + A _9_ + J_ + _1_) M 

K^/2- 3x j o o 2K ^| K 3 / 2 [ 0 0 SK^ 4K. 7 ^ 2 2K 2 Sx 


3x 


r A /I . l. r l/2v 

E 1 c 0 (- + ^ ) 


K 


3/2 


= . 0 , 


(3.19) 
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By adding and subtracting identical terms, thus adding a well chosen 
zero to (3.19), its form can be changed to 


r) 


1_ E ° C ° (2K 4K 3/2 ' . _ p , 1 + _L_x 3K 3_ 

3x „3/2 + o C o \ r 7/2 T OT ,2' 3x 3x 


K 


2K 


a ,1 , 1^1/ 2. 

E 1 C o ( k + ? 5 


K 


3/2 


= 0. (3.20) 


Using (3.17), we may substitute for E q in (3.20) and determine the 
integral of the middle term, obtaining 


fMo ln K - 


(3.21) 


By thus using (3.21), (3.20) can now be written as 


Thus 


3_ 

dx 


2^o c o ^ 


2K 1/2 - 3 
4K 1/2 + 2K 


2> + E 1 £ o C 


K 


5/2 + 2K } + 2^° C o ln K 


= 0, 


(3.22) 


3± a , 2K 1 ^ 2 » 3 
c (- 


3t A 


2 0 0 4K 1/2 + 2K 


+ E 1 C o^372 + 2K> + 2 E o c o K " c °* stant 


(3.23) 


follows from (3.22) . Again we choose to develop this constant at the 

A A 

locations where U — 0, K = K = 1, = E-^. Thus we obtain 


constant = §E Q £ o (- j) + + E 0 £ Q (o) 


or 


3 A A 1 A A. 

constant = - ^o c o^ 


(3.24) 


A . 

Notice the E n appearing in (3.24), something we have no expression for. It 

J- A 

is thus at present an unknown constant. To determine E^, we have to look at 
the equation of total energy for a capillary-gravity wave, given by 


E - E + E 

capillary gravity 



or 


1 2,2 1 2 1 2.2,.. g . 

-rpfa k + ^-pga = ypya k (1 + — p 

. * Z Yk 


Then, where U = 0, E = 


* * A 

E, k - k, a - a, so that 


But this is just 


£ 1 *2£2,_ g 

E « ypya k (1 + —p 

Yk 


A A A A 2 A 

E « E (1 + $ = E + e E- + e E 0 + 
o o 1 2 


Thus (3.27) results in the general condition that 


E ± « 0, i > 2. 


From (3.27), we have now that 


A A 

E- - E , 
1 o’ 


and (3.24) can be written as 


3 A 1 A A 5 * A 

constant = — {E c - fE c } = yS c ♦ 
2 o o 6 o o 4oo 


Combining (3.23) and (3.30), our result is -now 


2 E oV^l/2 + P + In K] + + 2R ) 


or 


\5/2 

E 1 ^ 

__ «, 

E. 


3K 5/2 m K - ^ l2 (- 2 4~-S ) 
4K 1 + 2K 

K 3/2 + 2 


Equation (3.32) is thus an equation for order e ^ in energy for z 
gravity wave meeting a current U. 


(3.25) 

(3.26) 

(3.27) 

(3.28) 

(3.29) 

(3.30) 

(3.31) 

(3.32) 
capillary- 
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To obtain the result to order e » we now write the expansion of the 
energy equation (3.3) in order thus 


3 r _ a , 3 1 , i _ * , 21 

{E o c o ( ,,„7/2 8K )} E o c o ( „„„9/2 


1 , 3| 


16K 


32K‘ 


16K 2 4K 5/2 2K 1/2 3x 


+ f- {E c 

9x 1 o 2K 


4K 


3 y , . _ * f 9 __1_ , 1 s 3K 

3/2 )} + E 1 c o ( 3/2 + 2 + i/ 2 > 9x 


8K' 


4K 2K 


+ {E 2 C o ( K + 2* ^ + E 2 C o ( ” 2k 2 ~ 4 k 1/2' 3x 


3 )|*- 0 . 


(3.33) 


Note that the E 2 terms have the same form as the order e case; the 
terras the same fora as the order t° case. Equations (3.16) and (3.32) 
specify what E Q and E^ are, thus a similar approach may be used here starting 
as before by multiplying through by thus making possible the combination 

of the last two terras in (3.33). This results finally in 


3*3/2 

2 s 


E o <k + 


J-291 _ r 3 5 2K 5/2 1 
96 " 16K 5/2 + 8K 4 
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8K 
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15 , 3 , 1 - .1 ~s2 

32K + 1/2 + 6 ln K + 2 * ln 


5-2 


1/3 
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ln 
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1+2 


1/3 


K 1/2 + 2 1/3 


1 - 2 1/3 + 2 2/3 
K - 2 1 ^ 3 K 1 ^ 2 +*2 2 ^ 3 


15*2 


1/3 


32 


VT 1 


arctan 
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2* /3 - 1 


- arctan 
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1/2 




l2 4/3 - K 1/2 


(3.34) 



Thus equations have been determined for Pq (3.16), (3.32), and E 2 (3.34), 

and hence we may write the energy to order as 

E - E q + eE 1 + e 2 E 2 , (3.35) 


knowing only the dimensionless parameter K. Thus rewriting (3.8) as 


or 

we can define U 
get 


T<" + #) - 1 
k C C 


k(¥ + •?■) = 1, 

c c 


(3.36) 


U __ 
c 


c 

■a-, and substitute into (3.36), rearrange, and 


v - — ± — 

■ U + c * 


(3.37) 


One possible solution scheme then, is to numerically solve for K for 
various values of U and c, and then use (3.16), (3.32), (3.34), and (3.35) 
to obtain the energy value for each particular value of U and c chosen. 
This is easily done by computer. Note in (3.37) that singularities exist 
at values where U + c approach zero. For these singularities, K values 
become infinite. 

Computer analysis can now easily be done and checked with laboratory 


results . 
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PART II 


EXPERIMENTAL STUDIES 

Xn this part of the report, the wave measurement system was first 
discussed. The system was used to study wave decay and capillary- 
gravity interactions. Some preliminary results of wind waves are also 
included. 
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1. INTRODUCTION 

This section deals with experiments on Gravity and Capillary waves 
relative to improving the understanding of how these wave motions are 
influences by surface currents. It Is known that active microwave back- 
scatter results primarily from the capillary waves via the mechanism of 
resonant or Bragg scattering. The relation between the ocean waves that 
are the primary contributors to the radar cross section and the incident 
radar frequency is shown in Figure 1.1. This figures was developed for a 
radar depression angle (NADIR) of 45°. It is evident that short wave- 
length or capillary waves contribute the most to the radar cross section. 
Moreover Huang (1972) has also shown that the shorter waves are more strongly 
influenced by currents. Therefore, there should be a strong correlation 
between radar cross section and ocean surface currents. 

The dependence of radar cross section on surface current can be expected 
from this simple argument. However, these has been considerable difficulty 
in making measurements of the short or capillary waves, and thus difficulty 
in relating the effect of currents via the change in capillary wave spectrum 
to radar cross section. Indeed, until recently and with the exception of 
microwave measurements of the ocean wave spectrum, the only method of making 
capillary wave measurements was with small wave height probes. McGoldrick 
(1971) has developed a capacitative wave height probe, specifically designed 
for capillary wave measurements. Although the system has been used for 
capillary waves, there has remained some question as to the validity of using 
a probe for the measurement of short wavelength waves. This is primarily 
because any probe will have a . meniscus on it due to surface tension and 
because capillary waves are also driven by surface tension. Therefore it was 
not clear that any effect due to the meniscus could be removed without 



altering the wave under consideration. Because of this, various optical 
systems have been proposed to measure capillary waves or to verify the 
accuracy of the wave height probe. Recently Tober, et.al.Cl973) and Sturm 
and Sorrell (1973) reported optical wave measurement systems. The data of 
Tober, et.al.(l973) reveals higher frequency components than measured with 
a wave height probe for the same conditions. Moreover, Sturm and Sorrell 
(1973) made a careful comparison of the response of their optical system 
with a capacitative wave probe of the kind developed by McGoldrick (1971). 
These results show a gradual roll off of the probe sensitivity with increas- 
ing wave frequency or decreasing wave length. Based on these results, the 
practical limit of the wave height probe is for wave lengths greater than 
2-3 cm. This is a longer wave length than much of the wave spectrum of 
concern here, and thus the use of an optical wave measurement system is 
required. . 

The system employed by Sturm and Sorrell (.1973) used a one dimensional 
position sensitive diode. While this diode is commercially available, it 
restricts the system to measurement of one dimensional wave trains. Most of 
the data of interest involves 2 dimensions or random wave fields, and thus a 
different diode array or matrix is needed. Part 2 of this report describes 
the 2 dimensional diode matrix that is used and the necessary electronics 
developed for data acquisition with this system. 

The following parts of the report deal with capillary wave experiments. 
An important question is how long the capillary waves live after ' generation. 
While this is well understood when no surface current or possibly when a 
constant current is present , there have been no data for variable currents 
where wave-current interaction is expected. In addition the presence of 
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any surface film can alter the capillary wave damping. Because of this an 
investigation of capillary wave decay on constant and variable currents as 
well as an investigation of surface films was undertaken. Part 3 reports 
the results of these experiments. 

It is well documented that microwave backs catter occurs by the mechanism 
of resonant or Bragg scattering from the capillary waves* The presence of 
the longer waves or the ocean swell is not negligible, however, and there has 
been considerable work including this in computations of radar cross sections . 
The more successfull attempts have been with the so-called "composite model’ 1 
where the swell is considered to simply tilt or alter the orientation of the 


surface containing the capillary waves. This model requires a knowledge of 
the wave height spectrum for the capillary waves and an important question is 


the influence of the ocean swell on the distribution and amplitude of the 


capillary' waves . In simplest form this is the question of any interaction 
of a long or gravity wave with a capillary wave. To investigate this an 
experimental study of the interaction of gravity and capillary waves was . 
undertaken. The results of these experiments are given in part 4 of this 


report . 
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Figure 1.1 Relation between radar frequency and scattering ocean waves. 
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2« Wave Measurement System 

2.1 Introduction 

The overriding objective of this work is to aid in the development of 
remote sensing techniques to the point where the wave characteristics, in 
particular wave energy or slope spectra, can be used to infer ocean currents. 
This requires work in two areas (1) development of remote sensing techniques 
to measure the ocean wave characteristics and (2) relating these wave 
characteristics to the mean current. Most of the work reported here is 
devoted to advancement of knowledge relative to problem (2), which is relating 
the wave characteristics or spectra to the mean current that is present. 
However this cannot be accomplished without some knowledge of remote sensing 
techniques. 


Most of the nr p sent 1 v c namn+rt — - j-i — j-. — , • 

* - - -j v ’— * - otiioing me i.iiuuti3 «mpj.oy tne reflection 


of some form of radiation from the ocean surface. This reflection appears to 
be related to ocean surface roughness or to the slope distribution of the 
ocean surface, rather than to the height distribution of the surface. 


However most of the presently employed oceanographic instruments measure wave 
height and not wave slope. ^ For a single wave there is a direct relation 
between wave height and wave slope, however for a random wave field there is no 
such relation, although there are empirical equations that are sometimes used. 

It is for this reason it is considered desirable to develop instrumentation to 
measure wave slope and to try to relate this measured wave slope spectra to mean 
water current. Such instrumentation has two uses in the remote sensing program. 
(1) It can be used to check or verify the data obtained by various remote 
sensing techniques. For this purpose the instrument should be able to measure 
ocean wave slope spectra. (2) It can be used to obtain laboratory data which 
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are used to relate measured wave slope spectra to known mean water currents. 
For this purpose the instrument must have a high accuracy. 


Since our program is primarily an investigation of the relation of wave 
characteristics to the current , the development of an instrument to make 
laboratory measurements has been our major purpose. However, when possible 
we have also kept the requirements of objective (1) in mind and tried to 
develop a system that can be modified for ocean measurements. 

It should also be noted that much of the ocean’s surface roughness is 


due to the capillary waves and that the longer gravity waves which contain 

most of the energy only weakly affect the surface roughness. Therefore the 

most useful information is a relation between capillary wave characteristics 

and water currents. Any wave measurement technique that requires the imersion 

a 

of a probe through the water surface will have /meniscus formed on the probe 


due to the surface tension of the water. Since this surface lcuSiou is a-iao 
the driving force for capillary waves there has been considerable discussion 


by many people as to the accuracy of any probe technique for measuring 
capillary waves. Recently Sturm and Sorrell <1973) have shown experimentally 
that the meniscus does cause significant error when measuring the wave height 
of short capillary waves with a probe. Therefore it is desirable that the 
measurement of wave slope be made without a probe being immersed through the 
surface. Finally, since the capillary waves have very short wave lengths, 
good spatial resolution is required. Thus the requirements for the wave 
measurement system are that it (1) measures more slope to a resonable accuracy 
(+‘ 2%) (2) that the spatial resolution be at least 0.5 mm and that the frequency 
response be at least 100 Hertz. (3) that no probe sticking through the surface 
be required. (4) If practical, a system that can be adopted to ocean surface 
measurements be employed. 
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The system developed to meet these requirements is an optical wave slope 
measuring system, and at present has met all criteria mentioned except (4). 
Moreover, at this time there do not appear to be an insurmountable problems 
in meeting criteria (4) . 

2.2 Principle of Operation 

Because of the inherent disadvantage with the probe system, there have 

been many optical techniques suggested for the measurement of wave slope. A 

brief review of some of the most applicable techniques is given with a 

that 

description of the present system. Early techniques/were considered by 
Hulburt (1934), 6ox and Munk (1954), and Schooley (1954), used the reflection 
of ight to infer ocean surface roughness. While these techniques or modifica- 
tions of them have been considered as possible remote sensing systems , they 
are not suited for accurate laboratory measurements. Moreover, they are not 
generally considered desirable for remote sensing because of the inherent 
difficulties in data reduction. 

The use of the refraction of light at the air-water interface is a more 
accurate technique. In this method a thin beam of light is usually directed 
vertically through the surface and the angle that the light beam is refracted 
from the vertical is measured. This system was first proposed by Cox (1958), 
and he measured the refraction angle by intensity variation. He experienced 
some errors due to intensity variation due to factors other than refraction 
angle and also due to water level changes. Later Prettyman (1969) used 
refraction and recorded the refraction angle with a high speed movie camera. 

this technique is accurate, data reduction is tedious for simple wave 
systems, and impractical for a random wave field. Recently Tober (1973) 
reported a refraction system in which an optical arrangement removed the effect 
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of changes in water elevation. However the refraction angle is measured by 
a variable density transmission filter, which converts the refraction angle 
to a change in intensity. Thus any phenomena which results in intensity changes 
(spray, water impurities, etc.) will cause errors. This system is accurate 
for water level changes up to 6 cm. 

The present system measures refraction angle directly and thus changes in 
light intensity do not affect the data. The major requirement is that the 
distance of the detector from the mean water level must be large compared to 
changes in mean water level. 

In order to illustrate this the principle of operation of the system is 


described. Fig. 2.1 illustrates an idealization of a one -dimensional wave as 
traversed by the incident light beam. 


m-L -I « _V-J- 1 ^ Vv«1/-\T.T -t-Vi Q su^ac* 

iflc ucam culci o » Oj- k. -i. j j. i 

which has the instantaneous slope d£/dx = tan0^ at the angle 6^ with respect 
to the surface normal. As a result of refraction the beam enters the air at 
angle 8^ with respect to the same normal. Angles 8 1} © 2 and indices of refrac- 
tion, n x and n 2 , are related- by Snell's law of refraction, n ± sin B ± = n 2 sin 8^ 
Since 8^ and 8^ are measured from the surface normal, the difference angle 
a = 0 2 - 8 X or simply the angular deviation from the vertical is measure experi- 
mentally. The difference angle a can be related to the angle through Snellfe 


law to give 


cotCGp = (n 1 /n 2 ) esc a - cot a. 


( 2 . 1 ) 


The surface slope d^/dx = tan9 1 is simply the reciprocal of equation 
(2.1). Hence, by a measurement of a one can obtain the instantaneous surface 
slope at the point traversed by the light beam. Equation (2.1) is valid for 
any wave slope as long as 0^^ is less than the critical angle 0 c where total 
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internal reflection occurs. In order to measure the difference angle a, one 

relates the offset distance d of the refracted beam from the vertical to the 

distance Z of the distrubed surface to the photodetector as simply a = tan^Cd/S.). 

When £ >> a, where a is the total vertical surface displacement, Z may be 

considered constant. The present system requires that Z be constant, and thus 

for the method in its present form to be used for finite amplitude waves, the 

detector must be sufficiently far away that Z » a. This is the only restriction 

in wave amplitude that is inherent in the system. Previous results using this 

technique which were reported by Sturm and Sorrell (1973) employed a United 

Detector model PIN-LSC-9 photodiode which has an active area of 2.5 mm by 225 mm. 

Because of the narrow width of the diode only one-dimensional waves could be 

measured. If there is any cross motion the beam refracts perpendicular to the 

diode axis and does not fall on the active part of the diode.. In addition the 

length of the diode puts an upper limit on the maximum slope -and amplitude 

that can be measured. The requirement that Z >> a determines a minimum 

distance, Z • between the diode and the mean water level. However for a fixed 

, min > 

vertical refraction angle, ct, and diode length, there is a maximum length, ^ max J 

for which the light remains on 'the active part of the diode. Thus £ max is 

and 

determined by the diode length /maximum slope to be measured and is determined 

by the requirement that Z >> a. Therefore there is a restriction on both the 
minimum and maximum distance of the diode from the mean water level. While this 
arrangement permits laboratory experiments under some useful conditions (see for 
example Sturm and Sorrell (1973)) it basically restricts its use to low amplitude 
one-dimensional wave trains. It is clear that modifications of the system are 
necessary if it is to be used in studies relative to remote sensing applications. 
The required modification are in the diode matrix used to measure the refraction 
angle. 
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2.3 The Diode Matrix 

A new diode' matrix was devised to greatly increase the active area of the 
diode system. A large area position sensitive diode is not presently available, 
and even if it were it would be prohibitively expensive. The present approach 
has been to use an array of diodes, in which each diode gives a signal that 
depends only on the fact that light is incident on the diode. The diodes are 
arranged in an array or matrix to give the size active area desired, and suitable 
electronics employed to give a signal whose voltage is proportional to the 
position of the diode upon which the light is incident. 

The diodes utilized are Vactec S-150-LB diodes which have an active area 

of 15cm X 6.5 cm. A schematic of the diode arrangement is shown in Fig. 2.2. 

Two diodes are placed end to end giving the matrix a width of 30 cm and 18 

diodes units are placed parallel to each to produce a length of 40 cm to 100 cm, 

depending on the diode spacing. In order to measure slope in 2 dimensions the 

refracted light beam is divided by a beam splitter and directed onto 2 diode 

matrices , one which measures position in the x-direction and one which measures 

position in the y-direction. This arrangement is also shown in Fig. 2.2., The 

accuracy of the measured slope depends only on the number of diodes employed. 

For 18 diodes the accuracy is — = 5.5% if all diodes are used. The spacing 

-KXXXKM between diodes can be changed to produce an active area that is as long 

as desired, and thus the system only requires that the array be the minimum 

distance, l . , from the mean water level, 
mm 

These diodes are purchased "off the shelf" from Vartec electronics and 
are quite reasonable ($1.00 - $2* 00/each). Thus the dxode 1 array presents no 
particular problem. However the associated electronics to convert the signal 
from the diode to a voltage proportional to its position in the diode array is 
quite involved. This has been a major development item for the present investi- 
gation and is thus described in detail. 
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2.4 Electronics, for the Diode Matrix 

The desired behavior of the diode matrix electronic circuit is to convert 

the signal output of a single diode to a voltage that is related to that diode's 

position in the matrix. The conversion of many similar outputs to a single 

coded signal is best handled with digital logic, therefore the approach has 

been to amplify the diode output to a voltage compatable with a digital logic 

converted into 

system. By using digital logic elements the signal can be easily/to a binary 
code indicating the diodes position. The resulting binary code can be converted 
to an analog signal, if desired, by a digital to analog (D/A) converter. The 
circuit for each step is now described in more detail. 

■ The basic circuit for the diode output is shown in Fig. 2.3, the diode 
amplifier circuit. The diode (indicated by PC) output goes to the first stage 
of a dual operational amplifier package, IC1. The first amplifier is used to 
amplify the current out of the diode. That is, this amplifier is. set up for 
current amplification. The amplified current goes to R2 which is of variable 
resistance to adjust the voltage input to the second stage of the dual system. 
The current output for a fixed incident light intensity is not exactly the same 
for each diode and R2 is used to compensate for this such that the voltage into 
the 2nd stage is the same for each of the 18 inputs. This voltage goes into one 
side of the 2nd amplifier, which is wired to operate as a comparator. If the 
voltage is above the reference voltage, R, the output rises until it is clamped 
at 4.0 volts by Dl, a Zener diode. If the output is less than R the voltage 
out of the 2nd stage is less than 0.5 volts. R3 is simply,. a current limiting 
resistor. These voltage levels are TTL compatible and thus TTL elements can 
be used for the rest of the circuit elements. 

A reference level is employed to allow the diode matrix to be used with 
high background light levels. That is each diode has some current output when 
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there is background light and there is an additional output from the diode when 
the refracted light is incident upon the diode. Thus it was considered desirable 
to have a reference signal level, such that there is no signal output until the 
refracted light is incident upon the diode. Because the background light level 
is highly variable, it was desirable to have the reference level controlled by 
the diodes themselves. Fig. 2.4 shows the reference voltage circuit which is 
driven by 2 separate diodes at the top and bottom of the diode matrix. This 
circuit follows the background light level and supplies a reference output for 
ICl^ the diode comparators. This circuit has variable feedback through R24 so 
the amplification of the background light can be controlled. The circuit is 
designed to have a very low output impedance to prevent switching of the 
comparators to alter the reference voltage level. This circuit appears to work 
quite well and has successfully removed the background from direct f lores cent 
light which oscillates at 120 Hertz. 

The output from IC1 has a slow rise time relative to that normal for 
conventional digital logic elements. When these logic elements experience a 
slowly rising voltage they exhibit instability and oscillate while the voltage 
is between the 2 logic levels. For this reason the output of IC1 is directed 
through, a current limiting resistor into a Schmitt trigger, IC2. The Schmitt 
trigger shapes the output into a fast rising signal. This signal is then run 
in to an 18 input- to 5 bit binary converter made up of NAND gates and NOR gates. 
This circuit, including Schmitt triggers is shown in Fig. 2.5, as the 18 input 
to 5 bit binary converter. The 5 bit binary output is indicated as outputs A 
through E. The choice of 18 inputs is arbitrary as up to 31 diode input channels 
can be employed with the 5 bit binary system (1 channel must remain at zero). 
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The 18 input to 5 bit binary converter works for any signal input from the 
.18 diodes, however anytime the signal is removed the output would be 5 bit 
binary zero. It is desirable to have the output remain at the last diode 
position or logic level until the next diode is turned on by the incident light 
beam. This is accomplished by ICS, which is a Quadruple Bistable Latch. The 
Latch is triggered by any input and holds the signal input until a signal from 
another diode is received. The circuit is shown in Fig. 2.6, the 5 bit binary 
to analog converter. ICS is a Hex Inverter which inverts the signal 4 times 
and is used only to add sufficient delay for the Latch to operate. The result 
is that any diode output is held by the Latch until an output from a different 
diode occurs. This produces a signal which does not drop to zero when the light 
beam is between adjacent diodes. 

The 5 bit binary output from the Latches goes directly into a D/A converter, 
IC7 , if an analog output is desired. If binary output is desired it can be 
taken. directly out of the Latch. IC8 amplifies the analog output to the 
desired voltage level. 

In most cases only an AG signal is desired. That is any DC offset because 
the unrefracted light beam is not incident on the center of the matrix is not 
wanted. Therefore a separate AC coupling amplifier is used. The circuit is 
shown in Fig. 2.7, the output amplifier. The arrangement provides a direct DC 
signal and also an AC coupled signal. The operational amplifier, IC12, is used 
to provide a high inpedance source so AC coupling down to 0.1 Hertz is obtained. 

The output from either the AC or DC jack is made of a series of small 

steps. This is shown in Fig. 2.9, which is a record of the signal out for a 

single wave train. The steped output, which is due to the position of the 

probably 

individual diodes, could /cause errors in -a- spectral analysis of the signal. 
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For this reason an analog smoothing circuit is also included in the output 
amplifier. This is also shown in the circuit (Fig. 2.7). It is simply another 
operational amplifier used to smooth the output. Fig. 2.9 shows the DC output 
and the smoothed output for a single wave. The input impedance R16 and feed- 
back capacitance C4 are variable to produce the desired amount of smoothing. 

For the' data given in Fig. 2.9 R16 was 100K and C4 = .002 

In addition to smoothing, in some cases it is desirable to Integrate the 
wave slope to obtain wave height • While there are some scaling problems 
associated with obtaining wave height in this manner, see for example Sturm and 
Sorrell (1973), many situations arise when it is desirable to be able to see 
wave height versus time, even with the height axis unsealed. This can be 
accomplished by changing the input inpedance R16 and feedback capacitance C4 
to values which make the smoothing amplifier and integrator. Fig. 2.10 3uu*»a 
the integrated wave slope measured by the optical system and the wave height 
as measured by a capacitance probe. In this application R16 is 39K and C4 is 
0.5 yf. 


All circuits were made on printed circuit board which was etched from 
photographic templates that are made in our Laboratory. The 18 input to 5 bit 
binary circuit with Latch and D/A conversion was the most difficult circuit 
to wire. A copy of the template used to make the printed circuit board for this 
part of the circuit is shown in Fig. 2.8. A list of all parts used in the 
circuit is given in Table 2.1. 


2.5 Summary 

An optical system to measure random wave slope has been developed and 
operated in the laboratory. The system has spatial resolution less than 0.5 mm, 
frequency response greater than 100 Hertz and does not distrub the surface in 
any way. Output from the system has been used to measure wave slope and wave 
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height for single wave trains, as reported by Sturm and Sorrell (1973). 

The system is presently being used to measure the interaction of capillary 
waves and gravity waves and to study the interaction of capillary waves with 
mean water currents. The only meaningful data for real wave— current inter- 
action is to measure the wave slope spectra for a random wave field with 
current. The system has been designed to do this specifically and such spectra 
are presently being obtained. 



Table 2.1 - Parts List' 


Integrated Circuits 


Part 

Number 

' Function 

IC1 

yA747 

Operational Amplifier 

IC2 

SN7413 

Schmitt Trigger 

103,109 

SN7430 

8 input NAND gate 

IC4,IC10 

SN7402 

NOR gate 

ICS 

SN7404 

Hex Inverter 

IC6 

SN7475 

Quadruple Bistable Latch 

IC7 

MC1406 

D/A Converter 

IC8 

yA741 

Operational Amplifier 

1012,1013,1014 



Active Components 



Q1 

40409 

NPN power transistor 

Q2 

40410 

PNP power transistor 

D1 

1N747A 

3.9 Volt Zener Diode 

PC 

S 150 LB 

Vactec Selenium photovoltaic cell 

Capacitors - All Capacitors 50.WVDC 

ceramic unless otherwise specified 

Capacitors 

Value 


Cl 

100 pf 


C2 

1 yf 


C3,C5 

0.1 yf 


C4 

Integration 

Smoothing 

0. 5 yf 
0.002 yf 


C6,C7 

10 yf 25 WVDC electrolytic 
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Table 2,1 (cont.) 


Resistors - All resistors 1/4 watt unless otherwise specified. 


Resistor 


R1 ,R3 ,R5 

IK 

R2 S R14,R23 ,R30 

10K Trimpot 

R4 

IK Trimpot 

R6 

4.7K 

R7 

6. 8 meg . 

R8 

2 meg. Trimpot 

R9,R17 

2.2 meg. 

R10 S R11,R19,R20 

10K 

R12 ,R13 ,R21,R22 

1.5K 

R15 a R25 ,R26 

100 n 

R16 

Integration 

39K 

Smoothing - 

100K 

R18 

100K 

R24 

IK Linear Taper Pot 

R27 ,R28 

IK - 1 watt 

R29 

270 ft - 1 watt 
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Figure 2.1 Optical path and refraction angles for a one-dimensional wave 
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(A) Wave Height - measured by a capacitative probe 



(B) Wave Height - by integrating wave slope 

Fig. 2.10 Wave height as measured by a wave height probe and by the 
optical system (vertical scales are not identical). 



3. Wave Decay Studies 
3 . 1 Int roduct i on 

The investigation of capillary wave decay or attentuation rates on 
mean water currents has been nearly completed. Special attention was 
given to two separate conditions (1) when a constant mean water current 
is present and theoretical analysis predicts no energy transfer between 
the wave motion and the current, f . e. , there is no wave current interaction 
and (2) when there is wave motion with a spatially varying current and wave 
current interaction is expected from theoretical considerations (Phillips, 
1966). Generally one expects an exponential decay of wave energy, and 
simple tests of wave decay confirm this. However, with capillary waves 
the presence of any film on the surface may contaminate the surface and 
although the decay is still exponential it may be greatly different than 
that predicted for clean surfaces * McQoldrick (1970) — - well as Davies 
and Vose (1965) have reported results on this effect. 

Because the effect of surface films can alter the results the present 
work includes extensive tests for the presence of surface films . This work 
indicates when these films can be expected, the effect of the surface film 
on wave damping, and how surface films can be avoided in laboratory 
experiments . 

Another phenomena occurs in laboratory wave experiments, which is the 
occurrence of cross waves or waves that propagate perpendicular to the 
direction of the longitudinal wave under study. These waves are experi- 
enced in many wave studies when no currents are present and Mahony 
et at . (1972) have reported results from extensive cross wave experiments. 
The present study has also experienced cross waves under most conditions. 
Moreover the appearance and subsequent behavior of these cross waves is 
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not that expected, observed or predicted by previous work for waves 
motion when no current is present. The preliminary test results indicate 
that (1) cross wave current interaction can occur even when there is no 
energy transfer from the current to the longitudinal wave; and (2) the 
presence of a current can generate cross waves as the longitudinal wave 
decays. Both of these results have important applications in physical 
oceanography. The latter is especially important, because it indicates 
that a current can break up a regular longitudinal wave train into random 
motion even when there is no energy transfer from the current to the 
longitudinal wave train. 

The importance of both film damping and cross waves in wave decay 
studies with application to physical oceanography has caused both to be 
investigated in detail. The results of these investigations are reported 
here. 

3.2 Studies of the Effect of Surface Filins 
Numerous investigators have studied the problem of wave damping by 
viscosity and/or surface films. A recent example of this work is con- 
tained in a paper by Davies and Vose (1965) who measured the damping of 
capillary waves on both clean surfaces and those with controlled amounts 
of contaminations. In the absence of a surface film (clean surface), 
Davies and V-ose observed that the damping of waves resulting from viscous 
energy dissipation in the bulk of the fluid could be accurately predicted 
by the relationship (one-dimensional) 

E = E^ exp [ (-4v k 2 /c ) x] 
o r g 

where v is the kinematic viscosity, k the wave-number, and c the wave 

& 


f 


(3.1) 
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group velocity. This result is well known end can be obtained directly 
from a first order (linear) wave analysis (e.g , j Phillips, 1969). The 
wave energy E, also calculated from the linear analysis, is given by 


E = , P q2a . 2 : (3.2) 

2k 

where p is the fluid density, a the wave amplitude (sinusoidal wave), 
and o the radian frequency. In deep water (d^ > X/2) o can be expressed 
in terms of k by the relationship a 2 = gk + yk 3 , where g is the 
acceleration of gravity and y the kinematic surface tension. With this 
equation (3.2) can be rewritten in a more useful form 


E = £<i*2!x ( i + _a_ > = £M>!x ( i + 3 ) 0.3) 

2 yk 2 2 

where £ = g/yk 2 is dimensionless. The term ak is the maximum slope of 
the wave £ = a sin (kx - at), i.e., ak = (dc/dx) max ; furthermore it 
follows that (ak) 2 /2 is equivalent to [ (dC/dx) ms )] 2 where rms is a 
root-mean square average of the wave slope. For short water waves 
(X < 1.7 cm) 6 is less than one and in the limit as X -+ 0, B 0. For 
short capillary waves (in these experiments) X < 1 cm and 3 < 0.345 
equation (3.3) can be approximated by 


* ; " <£> ! 

rms 


(3.4) 


allowing equation (3.1) to be written as 


2 

E = py (^“) exp [(-4v k 2 /c ) *3 j 
0 rms & 


(3.5) 


which implies 
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= ^ exp cc_2v k2/ y x] - 

rms ' ^rms 


(3.6) 


Thus, for short capillary waves, one would expect to observe a wave 
with initial slope (dc/dx)^ decay at a rate given by - , 


A 

v 


2v k 2 


g 


(3.7) 


which is the viscous attenuation coefficient. In view of this, a number 
of capillary wave decay studies were performed on tap water. Special 
care was taken to minimize the long term growth of a surface film (whose 
existence would invalidate the preceeding analysis). This was accomplished 
by a constant skimming process whereby water added to the bulk of the 
fluid was removed from the surface by a slightly submerged sharp-rimmed 
standpipe. The rate of surface runoff was maintained at approximately 
10 cc/min. 

In each instance the experimentally measured decay was exponential 
with a logarithmic decrement (attenuation coefficient) A^ as given by 
equation (3.7). The logarithmic decrement is defined as the natural 
logarithmic ratio of any two successive wave amplitudes (slopes) divided 
by the wave length A = £n(a^/a^. + ^)/X or. simply A = £n(2)/d2 where d^ is 
the distance required for the wave amplitude (slope) to decay to one-half 
the original amplitude (slope). The results of one of these decay trials 
are presented in Table 3.1 and Figure 3.1. 
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Table 3.1 Summary of wave conditions for decay 
data presented in Figure 3.1 


Measured Calculated 


T = 23° C 


R = 321.59 
0) 

y = 71.29 

3 , 2 

cm'Vsec 

A = 0.0383 

V 

X = 0.733 

cm 

d 2 = 18.1 cm 

t = 28.21 

ms 

A _/A =3.17 

f V 


c = 25.98 cm/s c = 26.91 cm/s 

A = 0.0378 cm" 1 c = 36.11 cm/s 

s 


Although the surface of the wave tank was constantly skimmed during 

all wave decay experiments, the possiblity of a surface film existed. 

A practical upper bound on film .damping , which has been confirmed by 

experiment (McGoldrick, 1970), is that due, to a rigid (inextensible or 

close packed) film acting in conjunction with viscous decay. The 

logarithmic decrement for this film is given by A f = (avk 2 /8 ) 1;/2 /c g 

or a factor 2" 5 ^ 2 larger than A y due to viscous dissipation alone. 

For all waves with R > 32, i. e. , capillary waves with X > 0.2 mm and 

w 

a < 10 5 , A _ is greater than A . For these experiments the Reynolds 
r v 

number R varies between 200 and 350 which implies 2.5 < A^/A^ < 3.5; 
thus the total decrement A T would be 3.5 A^ < < 4.5 A^. McGoldrick 

(1970) experimentally obtained for unskimmed tap wat^r a total decrement 
of A t = 5.77 a v for a wave Reynolds number R^ = 450, which is approxi- 
mately what one would obtain from rigid film attenuation combined with 

viscous dissipation at this Reynolds number. 



RMS WAVS SLOPE (d^/dx) x 3.93 


1-0 


.8 



.4 


.2 


.1 



Figure 3.1 The measured decay of capillary waves on a static surface 
(U = 0) (the solid line is the theoretical prediction) 
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These calculations demonstrate that film damping, if present, should 
be easily observed experimentally in the form of increased damping. The 
data, however, show no such effect since within the experimental error 
(less than + 10 percent) the measured decay' rate was accurately predicted 
by viscous dissipation alone. Therefore, in view of the precautions 
taken to maintain a clean surface, i . e. j constant skimming, and no support 
in the decay data for increased damping-, it was concluded that little or 
no surface film was present. 

3.3 Data Acquisition Techniques 

The data for these and all following experiments were obtained as 
follows. The wave slope was measured as described with our optical 
system. Changes in wave slope resulting from small variations in wave 
amplitude were minimized by averaging over several periods (n > 10) with 
a true rms voltmeter (DISA 55D35), wave maker period (t) was measured by 
a digital counter (General Radio 1911) and static surface tension (T 1 ) 
by a du Notiy tensiometer. Wavelength was measured by comparing the 
phase of the water wave, to that of the function generator (Wavetek 134) 
as viewed on a dual trace cathode ray oscilliscope (Tektronics 7704). 

The position of the incident light beam was set at a distance from the 
wave plunger where the phase of the measured water wave was judged by 
eye to be the same as that of the function generator. The light beam 
was then translated a total of five wavelengths (5X) and the average 
value obtained was chosen as the actual wavelength. This averaging 
technique was needed because of the difficulty in determining when 
exactly 2tt radians of phase, shift had occurred; in other words, if the 
error in determining the endpoints of the phase comparison is spread 
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over several wavelengths, a more accurate measure of X can be obtained. 

This averaging process was justified on the grounds that any small 
variation in X over five wavelengths, as a result of a possible change 
in phase velocity, would be insignificant compared to an error in phase 
measurement, even though the measure of one wavelength might appear more 
desirable. The experimental phase velocity c was simply calculated from 
c = X/t and A the logarithmic decrement, measured from the slope of the 
experimental decay data when plotted on a semi-log scale. 

The relationship between radian wave frequency a = 2nf and the radian 
wave-number k = 2 tt/t for waves m deep water can be expressed as 

cr 2 = gk + yk 3 (3.8) 

Since the phase velocity is related to the frequency and wave number by 
c = a/k > equation (3.8) gives the phase velocity as 

c 2 = g/k + yk ... . (3.9) 

In addition to the phase velocity c 9 it is .convenient at this point to 
compute the group c g = • Thus the group velocity is 

c g „ g . + * L± [gk + Yk 3 r l/2 (3.10) 

The calculated values of phase and group velocity (c and c^, respectively) 

were obtained from equations (3.8) and (3.9) using the static experimental 

value of surface tension: the wave Revnolds number R = (— ^r) was calculated 

w 

from the measured wavelength and frequency (a = 2tt/t) and the tabulated 
value for viscosity at the measured temperature. 
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3.4 Capillary Waves on Constant Currents 

Following these experiments a study of capillary wave dissipation 
on a constant current was undertaken. Modifications to the system 
required the addition of a recirculating pump and a section of the test 
channel with a reduced area 20 cm wide by 2 to 4 cm deep in order that 
current velocities as great as 30 cm/s could be obtained from moderate 
flow rates (Q <_ 2 L/s). The water current velocity was measured with a 
hot-film anemometer probe (developed by us) powered by a const ant- tempera- 
ture anemometer bridge (DISA 55D01). Within the resolution of the 
anemometer (+^ 0.5 cm/s), the velocity profile in the reduced area test 
section was uniform. The usual mechanical wave maker was replaced by a 
pulsed air wave maker which produced waves without physically contacting 
the water surface. This device consisted' simply of a narrow slit 2 mm 
wide and 20 cm long through which air pulses , generated by a large 
permanent -magnet loudspeaker, could be directed at the water surface. 

This change was necessary because a mechanical plunger stagnated the 
surface flow and generated a train of ripples, fixed on the oncoming 
stream, which interfered with the capillary waves under study. 

With the modifications described above, wave decay measurements 
with constant velocity currents directed against the direction of wave 
propagation were carried out. In a majority of these trials, however, 
increased cross-wave content, which became more severe as one moved away 
from the wave generator, was observed. In addition, the ratio of cross- 
wave slope to progressive wave slope at any fixed point remained nearly 
constant over a wide range of wave generator amplitudes . This feature 
is not observed in wave decay studies on a static tank where a neutral 
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stability- effect was observed. That* is, there exists a wave generator 
amplitude which would not spontaneously generate cross-waves but would 
sustain (and amplify) externally introduced cross mode distrubances . 

Below this critical (and frequency dependent) amplitude cross-waves 
cannot be maintained, while above it cross-waves are always present. 

The existence of these points of neutral stability has been observed by 
other investigators and several theoretical and experimental papers have 
been written about cross-waves . 

One of the most comprehensive studies to date is presented in a 
two-part paper by Mahony et at. (1972). . Mahony surveyed previous theories 
and proposed that there are two modes of cross-wave generation one near 
and another more distant from the wave generator, while Barnard and 
Pritchard conducted experiments with cross-waves generated by a mechanical 
(flap type) wave generator on a still (zero current) surface, these 
experiments document a number of characteristics of cross-waves, including 
neutral stabilities and growth and decay rates. From this experimentalist's 
viewpoint, probably the most disturbing feature of cross- waves as reported 
by these investigators is that cross -waves are not stationary in either 
time or space. Rather they grow and decay at slow rates with time constants 
on the order of 75 seconds when the wave generator is producing longitudinal 
waves with a fundamental frequency on the order of 5 Hz. 

When a mean current was present, no condition of neutral stability was 

v 

observed. It was first hypothesized that the pulsed air or pneumatic wave 
generator was the source of this increased cross-wave content. This 
hypothesis was tested experimentally by measuring the cross-wave content 
(in two separate trials) of two identical trains of capillary waves 
(t = 35.90 ms, X = 0.928 cm) on still water, one generated by the mechanical 
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plunger and the other by the pneumatic generator. The amplitude of each 
generator was adjusted to give the same rms slope at a distance of 5 cm 
from the plunger and the variation in rms slope of the longitudinal wave 
was measured in the transverse direction (parallel to the wavefront) as 
a- function of position from the wave generator. The results of these 
measurements, with the maximum transverse variation in wave slope 
expressed as a percent of the total wave slope, are presented m Figure 
3.2. Initially the mechanical generator exhibits a lower percentage of 
slope variation in the transverse direction; however, after about 10 cm 
the percentage of cross modes from either of the generators is essentially 
the same with an equilibrium value between 30 and 40 percent. It is 
concluded, therefore, that the cross mode content initiated by either of 
the two generators is comparable, and that the increased cross mode wave 
content observed when a mean current is present cannot simply be attributed • 
to the pneumatic wave generator. Rather it appears as if the current 
contributes to the growth of cross modes by some mechanism which is not 
presently understood. Figure 3.2 also presents a plot of the growth of 
cross modes when a current of -8 cm/s is present. The initial growth 
rate is nearly linear with distance with a final equilibrium value 
approximately double that observed on a static surface. 

3.5 Wave Decay Measurements on a Constant Current 

The apparent rate of wave decay is related to the cross mode content. 

When this content was low, the measured decay rate was that predicted by 

viscous dissipation alone. The decay was exponential with a logarithmic 

decrement A = 2v k 2 /c based on the actual (local) wave-number (cf. j 
v g 

Section 3.5) and group velocity. A typical plot of one such run is 



RMS WAVE SLOPE (d£/dx) x 3.93 



Figure 3.3 The measured decay of capillary waves on a current of -11.5 cm/s 
(cross-mode wave content is low) 
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presented in Figure 3.3 with accompanying Table 3.2. For the given 
conditions one concludes that film damping is not present by arguments 
previously outlined; i.e.j predicted film damping for R^ = 210 would 
be on the order of 2.5 times as great and the data do not support this 
fact. When cross mode content is high, the apparent decay rate is not 
exponential throughout. Rather, the decay rate is initially exponential 
in the vicinity of the wave generator, but about 20 to 30 wavelengths 
from the wave generator the decay data become erratic with large 
variations from the initial exponential decay rate (Figures 3.4, 3.5). 
However, from an observation of the data it appears that a lateral 
transfer of the progressive wave energy takes place; that is, an increase 
in energy at one traverse location (at a fixed distance from the wave 
source) is balanced by a decrease at another. 

Table 3,2 Summary of wave conditions for decay data 
presented in Figure 3.3 

Measured Calculated 


R = 210.6 
0) 

= 0.041 cm"* 

= 16.9 cm 
Af/ A v = 2.56 

= 24.641 cm/s 
U = - c = 12.37 cm/s 


T = 24. 5° C 

y = 71.3 cm 3 /s 2 

X = 0.986 cm 

x = 80.42 ms 

c = 12.26 cm/s 

A = 0.0461 cm" 1 
exp 

U = 11,5 cm/s 



RMS V7AVE SLOPE (di;/dx) * 3.93 



Figure 3.4 The measured decay of capillary waves along the tank centerline 
with a current of -8 cm/s (high cross -mode content) 




WAVE SLOPE (d£/dx) x 3.93 



DISTANCE FROM WAVE GENERATOR (CM) 

Figure 3.5 The measured decay of capillary waves one centimeter on either 
side of the tank centerline with a current of -8 cm/s (high 
cross-mode content) 
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RMS WAVE SLOPE (dc/dx) x 3.93 
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Consequently, spatial averaging of the decay data in the transverse 
direction removed the individual peaks and troughs and resulted in a 
smooth exponential decay. In practice this was accomplished by measuring 
the wave slope as a function of distance from the wave generator at five 
or seven transverse locations across the tank. (Usually this was done 
symmetrically about the centerline). The tabulated values for slope at 
each position x from the wave generator were then averaged (arithmetic) 
to obtain the spatial average which was taken as the representative value 
of slope, and hence energy, at this value x. The results of one such 
trial are presented in Figures 3.4 and 3,5, which illustrate the peaking 
described above, and the results of spatial averaging of these data are 
given in Figure 3.6 and in the accompanying Table 3.3. 

Table 3.3 Summary of wave conditions for decay 
data presented in Figure 3.6 

Measured Calculated 


T = 24° C 

y “ 71.4 cm 3 /s 2 

X = 0.026 cm 

t = 54. 59 ms 

c = 16.96 cm/s 

. A = 0.038 
exp 

U = 8 cm/s 


R = 264 
W 

A y = 0.0355 cm -1 

d^ “ 19.5 cm 

A ./A =2.87 

f v 

= 25.04 cm/s 
U=c^-c=8.08 cm/s 


As a result one must conclude from these experiments, including the 
special case where U = 0, the following: (1) in the absence of surface 



films the Kelvin dispersion equation (3.9) accurately relates phase 
velocity (or frequency) to wave-number; and, (2) the attenuation of 
waves occurs as a result of viscous dissipation as predicted by 
equation (3.6). In none of the cases, however, was energy transferred 
either to or from the current into the progressive (longitudinal) waves. 
However, Figure 3.2 indicates that energy is transferred from the current 
into the corss-waves or transverse modes. 

3.6 Capillary Waves on a Variable Current 

From theoretical considerations it has been suggested (Phillips, 

1966) that wave energy transfer to or from a mean current, t.e., wave- 
current interaction, should occur whenever waves propagate through regions 
where variable currents exist. This means that wave-current interactions 
depend on current gradients rather than simply the existence of a current. 
The previous experiments on constant currents support half of this argu- 
ment, namely, when no current gradient is present no wave-current inter- 
action occurs. To verify the second half of the argument, wave decay 
studies on a variable current were undertaken. In this section the 
results of the study are presented. 

The measurement of capillary wave decay on a variable current 
proceeded in a manner indentical to that of the previous constant current 
experiments with the exception of wavelength measurements. These were 
obtained indirectly as follows. From equation (3.8) one obtains 

°o “ (k o c 0 )2 = gk o + Yk o (3,1 

where k Q , c Q , and a Q are the wave-number, phase velocity, and frequency 
of waves on a zero current. 
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From this it follows that 

kg = (r + s) 1 ^ 3 + (s - r)^ 3 , 

where 

r 2 = (o 4 /T t2 )/4 + (g/T * ) 3 /27 , 

and 

s = o 2 /2T' 


(3.12) 


i 

» 

{ 

i 


The dispersion equation (3.9) can be used with an equation requiring 
the conservation of waves (e.g. Phillips, 1966) which for steady waves 
on deep water reduces to 

Oq = constant = ^o C 0 ” = k(c^ + (3.13) 

These equations (3.9) and (3.13) are then solved simultaneously for wave 
number resulting in a cubic equation of the form . 


(k Q /k) 3 + a 1 (k Q /k) 2 + a 2 (k Q /k) + a = 0 , 


where the coefficients, a^, a 2 , and depend only on the known quantities 
U, Cq , k Q , y and g are given by 

2Uc n + g/k. 


a l = - ( - 


L > ’ a 2 = (U/c 0 ) 2 , a 3 = yk 0 /c 2 . 


This equation is solved explicitly for k resulting in an expression of 
the form 


k = k (k Q , c Q » U, y, g) . (3.14) 

Thus, for each value of current U, one can obtain a corresponding value 
for k and hence for the wavelength X. 

Figure 3.7 is an experimentally determined plot of the current used 
in the following wave decay study. By convention, the current is 



CURRENT VELOCITY, U (CM/S) 
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RMS WAVE SLOPE (d£/dx) ><3.93 


1.0 



Figure 3.8 Decay of capillary waves on a variable current 

(of . , Figure 3.7) (spatial averaging was required) 


r 



considered negative since it opposes the direction of wave propagation; 
also, the magnitude decreases linearly as one moves away from the wave 
generator (located at x = 0). 

Waves with a period of 92*99 ms (10.73 Hz) were generated on the 
-17 cm/s current; the waves propagated through the variable current region 
and were finally dissipated in the -5 cm/s constant current region. The 
measured decay of these waves is presented in Figure 3.8. 

The data presented are the spatial average of seven traverses of 
the test section as cross wave content was- high due in part to the fast 
current and relatively short wave period. 

The dashed line on Figure 3.8 represents the expected decay from 
viscous dissipation alone. This was obtained as follows. At each data 
point the wave-number k was calculated from equation (3.14) and the 
measured current velocity. This local wave-number was then used to 
compute the viscous decay decrement at each point; the dashed line is 
the tangent" of these local decay rates. It can be seen- from the figure 
that the initial and final decay rates (slope of the dashed line) are 
tangent to the data. ' This is consistent with the previous findings for 
wave decay on constant currents since at either end of the variable 
current section constant current conditions exist. However, the over-all 
predicted rate of decay is too slow to explain the decay observed in the 
variable current region. 

An analysis which includes the effect of a variable current on the 
wave energy is given by Phillips (1969). For a one-dimensional train of 
waves propagating in the x direction, he shows that the conservation of 
the fluctuating component of the wave energy can be expressed as 

X5/ 
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4— [E (0 + c ) ] + S 42- = -e , 

dx g xx dx 


(3.15) 


where E is the wave energy, U the current velocity, c the wave group 

velocity, S the radiation stress in the x direction, and e the rate of 

energy dissipation per unit volume due to viscosity. From a linear 

analysis, he also shows that e = 2yyk 2 (ak) 2 - 4vk 2 E and 

S = E ( ~ + - r ) where 8 = g/yk 2 (introduced in equation (3.3)) is 

xx 2(1 +3) 

a dimensionless parameter. For this experiment 3 reaches a maximum value 
of 0.98 which implies that the acceleration of gravity is no longer 
negligible and must be included in the analysis. As a consequence the 
radiation stress S and the group velocity c are expanded in a power 

XX 2 

series with 6 as an ordering parameter. Neglecting all but the linear 
terms, one obtains first order perturbation relations for S xx and c g 
given by 


and 


S = E (3/2 - 8) 
xx 


c = (3 - 8/2) c/2, 
g 


(3.16a) 


(3.16b) 


Thus equation (3.15) can be written as 

4- CE(U + (3/2 - 8/4)c)] + [(3/2 - 3) + 4vk 2 ] E = 0 . (3.17) 

dx 

Equation (3.17) can be put in the form 


f + p(x) E = 0 , 


(3.18) 


where 


p(x) = [(5/2 - 6) ^ + (3/2 


6/4)' 4 s - + 4vk 2 ]/[U + (3/2 
dx 


- B/4)c] . (3.19) 
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Equation (3.18) has the general solution 

E = Eq exp (- | p(w) dm , - (3.20) 

0 

where E^ is the wave energy at x = 0 and p(x) represents a generalized 

decay coefficient for capillary waves on a current. It can be readily 

seen from equation (3.19) that when the current is constant dU/dx and 

dc/dx are zero, and p(x) reduces to the constant 4vk 2 /(U + (3/2 - 6/4) c). 

Since capillary wave energy is proportional to rms slope squared, the 

decay decrement for wave slope is half the above or A where 

U + (3/2 - 6/4) c = c , the group velocity. 

S 

For a variable current (3.19) is a function of U(x) which must be 
determined experimentally. In view of the difficulty in obtaining a 
compact functional form for equation (3.19) and the approximations 
implicit in equation (3.17), no closed form integration of equation 
(3.20) was attempted. Rather, two graphical techniques were employed. 

In the first case p(x) was plotted on a Cartesian scale", and the integral 
of p ( x ) was simply obtained by counting squares under the curve. In this 
manner a table of values for the definite integral J Q p(w)dui was obtained. 
Subsequently, E/E Q from equation (3.20) was evaluated at each data point 
and plotted on a semi-logarithmic scale. 

Within a constant, the results of these calculations are presented 
as the solid line in Figure 3.8. Alternately, p(x) was evaluated at the 
values of x corresponding to each data point. The value of p(x) thus 
obtained corresponded to the energy logarithmic decrement or simply the 
slope of E(x)/E q at each data point. Thus, the solid line in Figure 3.8 
can also be interpreted as the- tangent to all the predicted decay rates 
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obtained from p(x). Table 3.4 summarizes the calculations involved to 
obtain p(x) given the experimentally measured values for the current. 


Table 3.4 Numerical calculations involved in the evaluation of 
equation (3.19) 


^cm/ s 

c. : cm/s 

l 

g/yk 2 = 6 

4vk 2 

dU/dx 

dc/dx 

p(x) 

2£n2/p(x) 

-17 

25.879 

0,238 

2.18 

0.1 

-0.9038 

0.1119 

12.379 

-16 

25.425 

0.269 

1.93 

0.666 . 

-0.25 

0.1496 

9.260 

-14 

24.625 

0.341 

1.52 

0.666 

-0.236 

0.1259 

11.00 

-12 

23.90 

0.434 

1.197 

0.666 

-0.19 

0.108 

12.829 

-10 

23.496 

0. 550 

0.943 

0.666 

-0.117 

0.0946 

14.648 

- 8 

23.175 

0.696 

0.746 

0.666 

-0.10 

0.0798 

17.38 

- 6 

23.012 

0.875 

0. 594 

0.666 

-0.011 

0.0706 

19 . 61 

- 5 

22,987 

0.977 . 

0.531 

0.125 

+0.0357 

0.0665 

20.86 

- 5 - 

22.987 

• 0.97-7- 

- -0.531 

0.0 

- - 0 .-0 

— 0-.-0222-. 

- 62.30- - 

From the data in Figure 

3.8 it 

can be 

seen that 

the decay rate 


throughout the variable current region is nearly constant, although the 
wave-number (k) and the current velocity (U) decreases rather rapidly. 
This effect can be explained by wave-current interaction theory through 
an examination of equation (3.19). It can be seen that the sum of three 
terms (5/2 - 8) , (3/2 - 3/4) , and 4vk 2 , jointly divided by the 

local group velocity, controls the rate of wave decay. At the beginning 
of the variable current region where U = -16 cm/s, the ratio of 
attenuation resulting from viscosity to that arising from the gradient 
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terms, i,e, j = 4vk 2 /[(5/2 - 0) + (3/2 - 0/4) ] is 

1.93 

R = —r? = 1.71, so that for the shorter wavelengths encountered the 

rate of decay is controlled by viscous dissipation. As a current 

decreases, however, the wave -number decreases such that at U = -12 cm/s, 

1.20 

R v g = = which implies that the gradient terms account for 

nearly half of the apparent rate of decay. For u = -6 cm/s, the ratio 

R = ■■ ■ = 0.558 and the effect of viscosity is overshadowed by the 

vg 1.066 J 

gradients. The net effect is that the sum of the three terms (5/2 - 0) 

d.c o 

(3/2 - 0/4) — , and 4vk z remain roughly constant, and hence the rate of 
wave decay also remains constant. The data in Figure 3.8 support this 
observation. 

It is important to note, however, that the increased rate of decay 
from the combined effects of wave-current interaction and viscous dissipa- 
tion is actually apparent rather than real. That is, the action of the 
radiation stress is to transfer (conservatively) energy from the wave into 
the current rather than to dissipate it as is the case with viscosity. 

This can be shown for the case of pure capillary waves, i.e.j 8 - 0, by 
considering equation (3.17) without the term (4vk 2 ) resulting from viscous 
dissipation, namely, 

^ [E (0 + 3/2 o) ] + | E H = 0. (3.21) 

For pure capillary waves, equation (2.12) reduces to 

c 2 = Y k 2 • (3.22) 


Combining equation (3.22) with (3.10), one obtains 
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c 3 U c 2 

£_ + iL- £__i = o 

c c 2 
0 0 0 


from which it follows that 


~ = -2c (U + 3/2 c) 

dx ax 


(3.23a) 


(3.23b) 


Multiplying (3.21) by c 3 and combining with (3.23b), one obtains 


dc 


3 ~- [E(U + 3/2 c)] - 3c 2 E(U + 3/2 c) -r” = 0 


dx 


dx 


(3.24a) 


or 


—■ [E(U + 3/2 c)/c 3 J = 0 


(3.24b) 


Thus, without viscous dissipation 


E(U + 3/2 c) , 3 Eo 

— = const = — — 

c 3 2 c 2 


(3.25) 


where subscript (0) indicates conditions when U = 0. 

For a single train of capillary waves, E = “ py(ak) 2 (from equation 
(3.3)). Substituting into equation (3.25) results in 


(ak) 2 

< a oV 2 


<— > 3 

V C 0 ' 


E_ 

-t~ ' E o 


2 U 

o + 3 ‘o 


(3.26) 


E/E~ as a function of U/c_ is plotted in Figure 3.9. It can be seen that 
0 u 

the wave energy E decreases as -U/c 0 tends towards zefo (U/Cq = 0). There- 
fore, the energy is not dissipated as a result of the wave-current inter- 
action. Rather it is simply transferred from the waves into the current. 
Thus, the vertical distance between the solid and dashed curves in 

..irw^r 
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Figure 3.8 is a measure of the total energy transferred from the waves 
into the current. Ideally, if the starting current conditions could be 
re-established, this portion of the energy would be returned to the waves 
It is concluded that this detail investigation of wave-current 
interaction with variable currents confirms recent analytical predictions 
that is the energy transfered from the wave to the current is exactly 
that predicted. However, it is probably more important to note that in 
a variable current the cross waves interact with the current to gain 
energy while the longitudinal waves interact with the current to loss 
energy. In physical oceanography one is concerned with the eventual 
wave state after propagation through, a variable current region and these 
experiments that after a short period in many cases the cross waves will 
be almost as large as the longitudinal waves. This is considered to be 
quite an important observation. 
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4. Experimental Studies of Capillary-Gravity Wave Interaction 
4.1 Introduction 

In almost all remote sensing techniques a large part of the 'data (usually 
back scatter at some angle) occurs or is caused by the short wavelength waves on 
the ocean surface. Therefore attention must be given the shorter waves and how 
they effect the data that is obtained. Since these waves have much shorter wave- 
lengths than the longer gravity waves, they may not be uniformly distributed on 
the longer waves, but may collect at the peak, trough, front or back of the gravity 
wave. In addition since the waves receive energy and decay at different rates, 
their relative amplitudes may change as the wave system propagates. Both of the 
factors, non-uniform distribution of the capillary waves and relative amplitude of 
the waves will alter the signal obtained by the sensing system. In the ocean there 
is, of course, a spectrum of wavelengths from the long gravity waves to the shorter 
capillary waves and thus we are considering the effect the gravity waves will have 
on the capillary waves. Any changes in the short wave distribution or amplitudes 
will alter a return- signal from these waves . * . .... _ _. 

In order to Investigate this a simplified experiment was devised and undertaken. 

This was to generate a train of gravity waves at a single frequency or wavelength 
and superimpose on these long waves a uniform series of much shorter wavelength 
capillary waves (also at a single higher frequency). The composite wave system was 
measured near the plunger and then measured at various positions of increasing 
distance from the plunger in an effort to determine how the capillary waves inter- 
acted with the gravity waves. This data was analyzed to determine any non-uniformities^ 
in the capillary wave distribution and to see if any energy was being transfered 
between the wave systems. The experimental results are plots of composite wave height 
at various distances away from the generation source. 



4,2 Experimental Technique 

The waves were generated by a single plunger that was driven by a signal 
which was the addition of the two desired waves. Two sine wave generators were 
used, their signals added electronically and amplified to drive a large permanent 
magnet loudspeaker. The Speaker was connected to a plunger which was used to 
generate the waves. Figures 4.1a and 4.1b show a typical signal feed into the 
amplifier to drive the plunger. This particular signal has a lower frequency of 
4 Hertz for the long wave and a frequency of 18 Hertz for the short wave. This 
signal record shown in Figure ^.JJb is .51 sec. m duration and the longer scale 
Figure 4.1a is 1.25 sec. in duration. Note that the amplitude of the capillary 
waves is considerably larger than that, for the gravity waves, the reason for this 
will be discussed later. This is the signal used to generate the data shown in 
Figures 4.3a to 4.3i. Figure 4.2 illustrates the speaker and plunger arrangement 
used t;o generate the waves. In order to assure that the plunger accurately 
reproduces the input signal to the amplifier a positive feedback circuit was 
employed. Figure 4.2 also shows this feedback circuit. 

After generation the waves were measured by .the optical wave slope system 
that we have developed, Sturm and Sorrell (1973). This system used the new 
electronic detection system and diode matrix that has been recently developed and 
was described in section two (2) of this report. This system measures wave slope 
directly and the signal must be integrated to provide wave height. The integrator 
used to accomplish this is also described in section two (2) and 
a discussion of the complete method may be found in the paper by Sturm and Sorrell 
(1973). One consequence of this procedure is that the amplitude scale is not the 
same for both waves, but varies directly with the phase speed of the wave. This 
real error however, as the phase speeds are known and if the absolute 



161 


ratio of amplitudes is desired, they r can be easily scaled. Moreover for the 
present experiments the amplitude scale was varied in order to show the waves as 
clearly as possible. Therefore, for all of the data presented, only the relative 
amplitude of the wave is presented and no absolute amplitudes are given. 

Both waves are thus generated at the plunger and are observed at various 
distances away from the plunger. The waves require some distance from the plunger 
to reach equilibrium before viscous decay starts. This distance is greater for 
the gravity waves than the capillary waves, and the intermediate position of 7 cm 
from the plunger is the first observation point in all cases. By the time the 
capillary waves reach this distance there has been some decay, however the gravity 
waves have experienced no amplitude loss. This is the reason the plunger motion 
for the capillary. waves (Figures 4.1a and 4.1b) is larger than that for the longer 
waves . 

A complete discussion of the relevant theory is present in the previous section 

of the report. However some of the experimentally significant results of the theory 

are given here. Basically the theory predicts a stronger interaction,- that is, the 

gravity wave will have a greater effect on the distribution and amplitude of the 

capillary waves, when the phase speed of the gravity wave is either equal to the 

3 

phase speed of the capillary wave or when . The greatest effect is 

3 

expected when because the energy transfer is at the group velocity and 

3 

the condition C 1 = — corresponds to equal group velocities. Table 4.1 gives the 

wavelength, frequency and phase speed for a range of capillary, gravity and 

capillary-gravity waves and is used to facilitate comparison of the wave systems. 

It is evident from this table that very high frequency capillary waves will be 

3 

required to meet the condition . The consequences of this will be 


considered later. 
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Table 4,1 

Wavelength, Frequency and Phase Speed in the Range 
of Waves Considered 


Wavelength 

Frequency 

Phase Speed 

X (cm) 

f (Hertz) 

c (cm/sec) 

30 

2.28 

68.6 

20 

2.80 

56.1 

15 

3.24 

48.7 

10 

4.01 

40.1 

8 

4.51 

36.1 

5 

5.90 

29.5 

4 

6.80 

27.2 

3 

. 8.30 

24.9 

2 

11.6 

23.2 

1.8 

12.8 

23.07 

1.6 

14.4 

23.07 

1.4 

16.6 

23.26 

1.2 

19.8 

23.74 

1 

24.7 

24.65 

.8 

32.8 

26.26 

.6 

48.4 

29.1 

.4 

86.3 

34.5 

.3 

131.3 

39.4 
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2.3 Experimental Results 

As previously mentioned the experimental results are plots of wave height 
verses time at increasing distances from the plunger. All data starts at an 


initial position 7 cm 

from the 

plunger. Data from 

3 separate conditions is j 

These conditions are: 




Run #1 

Gravity 

or long wave 

Capillary or short wave 


h = 

10 cm 

= 1.6 cm 


f i = 

4 Hertz 

f 2 = 18 Hertz 


c i = 

40 cm/sec 

C 2 = 23.9 cm/sec 

Run H 2 

Gravity 

or long wave 

Capillary or short wave 


x i = 

5. 1 cm 

X 2 = 1.1 cm 


f i = 

6 Hertz 

f 2 = 24 Hertz 


c i = 

30 cm/sec 

C 2 = 24.6 cm/sec 

Run #3 

Gravity 

or long wave 

Capillary or short wave 


x i = 

5.1 cm 

X 2 = 0.86 cm 


... f i = 

6 Hertz 

f 2 = 30 Hertz 


c i = 

30 cm/sec 

C 2 = 25.7 cm/sec 


The data for Run #1 is given in Figures **.3a to 4.3i, for Run #2 in Figures 4.4a 
to 4.4e and for Run #3 in Figures 4.6a to 4.6c. A comparison of the present optical 
wave measurement system with a conventional wave height probe can be obtained from 
Figures 4.5a and 4.5b, as these data (Figures 4.5a and 4.5b) were taken with a wave 
height probe. The results are now discussed in more detail. 

Run $1 

The data were obtained at the following positions from the wave source (plunger). 


4. 3a 
4.3b 
4.3c 
4.3d 


7 cm 
12 cm 
17 cm 
22 cm 


4, 3e 
4. 3f 
4 . 3g 
4.3h 
4 . 3i 


27 cm 
32 cm 
37 cm 
42 cm 


47 cm 
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The wave system data in 4.3a, 4.3b, and 4.3c- shows an even distribution and lends 
confidence that the waves were generated as expected. In Figure 4 , 3e there is 
some evidence that the waves are grouped more in the long wave trough, however 
4 . 3 f shows the waves grouped more on the long wave peak. Figure 4.3g and 4.3h 
show a similar result. At a position 47 cm from the plunger (Figure 4«3i) the 
capillary waves have almost completely decayed away. The capillary or short waves 
decay much faster with position than the long waves because they have traveled 
many more wave lengths in distance. This is another reason the capillary waves 
were generated with larger relative amplitudes. 

The apparent shift in position of the capillary waves from the gravity wave 

trough to peak and back is believed to be caused by the difference in phase speed 

of the two wave systems. By the time the group of capillary waves propagate to the 

next position the gravity wave has propagated approximately 3.4 cm farther than 

the capillary wave. This places the small waves at a dirrerent position on the 
wave 

gravity/. Indeed this illustrates that unless the gravity wave is sufficiently long 
for the capillary waves to die out between crests, a group of short waves, will 
change position on the long wave unless the phase speeds are equal. The present 
wave tank permits wavelengths up to about 12 cm before tank size becomes a problem. 
Thus the phase speeds must be equal or else a grouping of short waves will shift 
relative position on the long wave. 

Run #2 

The figures show data obtained at the following positions from the wave 
generation source. 

4.4a 7 cm 4 . 43 17 C m 

4.4b 12 cm 4.4e 22 cm 

4.4c 17 cm 
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Both the long wave and the short wave are at higher frequencies in an attempt to 
generate waves with the phase speeds more nearly equal. The main feature here is 
the rapid decay of the capillary waves , and by the time they have, traveled 17 cm 
from the plunger they have nearly vanished. Figure 4.4d is with an expanded time 
scale to verify this. Data given in 4.4e, which is 22cm from the wave source, 
contains virtually no short waves. 

Because detection of the short waves is marginal at a position 15.5 cm from 
the source, (Figures 4.4c and 4.4d), this data taken by the optical system was 
compared with that obtained from a conventional wave height probe. Figure 4.5a 
and Figure 4.5b (expanded time scale) shows data at identical conditions of Fain #2 ■ 
15,5 cm from the wave source. These data were taken with a capacitive wave height 
probe which had a probe diameter of 0.25 mm. In these data the capillary wave is 
completely obscured by the probe. This shows the high resolution of the optical 
system and the inherent difficulties associated with using a probe system to measure 
capillary waves. 

Run #3 

This condition has an even higher capillary wave phase speed in an attempt to 

determine if the waves group on. the long wave peak or trough. Data are shown for 

three positions: 

4.6a 7 cm 

4.7b 12 cm 

4. 6c 17 cm 

Figure 4.6a demonstrates that the waves are generated as desired and that there is 
ample capillary wave content. Figure 4.7b shows considerable capillary wave decay, 
and at a distance 17 cm from the wave source the capillary waves are no longer 
observable. This demonstrates the inherent problem with attempting to generate both 
a long wave and a short wave with the same phase speed. As the short wave length is 



(a) Trace duration = 1.28 sec 
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(b) Trace duration = 5.1 sec 
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Figure 4.1 Input to \7ave generator Run # 2 , = 4 S = 
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(a) 7 cm from wave source 



(b) 12 cm from wave source 


i 



(c) 17 cm from wave source 


Figure 4.3 Data from Run #1 






(g ) 37 cm from Vi 



(h) 42 cm from 







(a) Trace duration = 1.28 sec 



(b) Trace duration = 5.1 sec 


Figure 4.5 Data from wave height probe Run #2, 15.5 cm from source 
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(a) 7 cm from wave source 



(’e) 17 cm from wave source 


Figure 4.6 Data for Run #3 
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reduced in order to increase the phase speed, the wave decays very rapidly. If 
an attempt is made to produce shorter long waves with lower' phase speeds, the 
distinction between long and short waves is lost. 

4,4 Summary and Conclusions 

A simple experiment was undertaken to investigate the interaction between 
long gravity waves and much shorter capillary waves . The experiment was to 
generate both wave systems and to observe the effect of the long wave on the shorter 
capillary wave. Run #1 shows the difficulty in analyzing the results when phase 
speeds of the two wave systems are different. Run #2 indicates the rapid decay of 
the shorter waves with increasing frequency and compares the present optical system 
for data acquisition with the conventionsl wave height probe. This illustrates the 
necessity of using the optical system. Run #3 verifies the rapid decay of high 

frequency capillary waves and shows that it is probably not possible to generate a 

distinct long wave and short wave with the same phase speed and have the composite 
wave system last for any appreciable length of time or distance of propagation. 

For these observations it appears that a system of wind driven waves , or 

possibly plunger and wind driven waves are necessary for the experiments. If the 
short waves are wind generated, then it should be possible to continuously add energy 
to the short waves so they do not decay or at least do not decay so rapidly. With 

this approach the long waves could be plunger generated and the short wave generated 

and sustained by the wind. Under these conditions it should be possible to generate 
and measure long waves and short waves which have the same phase speed. 
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5. Preliminary Wind-Wave Measurements 

5.1 Introduction 

The previous sections present the results of measurements of waves generated 
by a wave maker or plunger. Waves generated in this manner have definite periods 
or wave lengths , and thus it is meaningful to describe the wave by these charac- 
teristics. However, wind generated waves are random and thus the only meaningful 
data are some form of wave statistics or average properties of the waves. Because 
radar backs catter is primarily influenced by surface roughness and also because 
-the high frequency components are important, results in the form of wave slope 
spectra are the most useful. These data can be obtained with the optical wave 
measurement system described in Part 2 of this section. Preliminary results of 
wave slope spectra for different conditions of current and wind speed are given 
in Figures 5.1. 5.2, 5.3. and 5.4. These results demonstrate a strong relation 
between current and slope spectra. 

5.2 * Data 

The data for the accompanying figures was obtained in the wind-wave tank at 
N» C. State University. This tank is 2 ft. x 3 ft. m cross section and has a 
fan capable of generating winds to 40 m/sec. In addition a pump system has been 
installed, to produce mean water currents. This system permits the generation of 
wind waves and the investigation of changes in wind-wave statistics due to a 

current. - 

The data were obtained by recording the output of the optical wave measure- 
ment system on magnetic tape. The tape system has a FM record and playback mode 
so that low frequency waves can be recorded. Data are recorded for a length of 
time sufficient to obtain the desired amount and then played back through an 
audio spectrum analyzer. Because the analyzer requires frequencies above 
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20 Hertz (audio frequency) the tape speed is increased by a factor of 16 when 
played back into the spectrum analyzer. This gives a low frequency limit of 
1,25 Hertz in the slope spectra. The output of the spectrum analyzer is then 
recorded on an x-y plotter. The data given in 5.1, 5.2, 5.3, and 5.4 are the 
analyzer output as plotted. These spectra shows considerable fine scale struc- 
true. In order to get a better picture of this structure several spectra are 
plotted on top of each other. This shows the ■ repeatibility of the spectra and 
also shows if a particular fine structure characteristic of the spectra or 
random. These slope spectra are basically plots. of the relative magnitude of 
the RMS slope within specified frequency bands. They are taken with a wind 
speed of 2 m/sec at zero water current and with a mean water current of 10 cm/sec 
in the direction of the wind. Figures 5.1 and 5.2 are wave slope spectra down 
the channel, that is, slope spectra in the direction of the wind or current. 

The shift of the spectrum is in quantitative agreement with the theory developed 
by Huang, et. al. (1972). Generally it shows that a positive current shifts the 
wave energy toward the -lower frequencies.-- The data given -in Figures 5.3 and 5.4 
are slope spectra across the channel. These have less energy content and, in 
addition, both the downwind. and cross wind spectra with the positive current show 
less energy transfer from the wind than when there is no current. Therefore it 
is not meaningful to compare the relative total energy content of the two spectra. 
The data show the expected strong dependence of wave spectra on current. 




Figure 5.1 Downwind slope spectra 

Wind = 2 m/sec Current = 0 cm/sec 



Figure 5.2 Downwind slope spectra 

Wind = 2 m/sec Current = 10 cm/sec 
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Figure 5.3 Crosswind slope spectra 

Wind = 2 m/sec Current = 0 



Figure 5.4 Crosswind slope spectra 

Wind = 2 m/sec Current = 10 cm/sec 
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